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ABSTRACT 

This thesis investigates stable, asymptotically stable 
and unstable dynamical systems from a topological point of 
view with direct application and interpretation to systems 
of differential equations. Knowledge of topology is not 
a prerequisite. Specifically, such concepts as Poisson 
and Liapunov stability as well as parallelizable and 
dispersive systems are investigated. 
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INTRODUCTION 



As it is currently available, stability theory of 
dynamical systems requires an extensive background in higher 
mathematics. The purpose of this thesis then is to make 
this theory accessible to the reader who does not have this 
background. It is assumed that the reader has taken courses 
in elementary differential equations, intermediate mathe- 
matical analysis, and linear algebra. All the topological 
properties needed are described in Section 1. 

Sections 1-3 deal with developing basic concepts of 
dynamical systems. Sections ^ - 6 use these concepts in 
studying various kinds of stable or unstable systems. Topics 
of recursiveness and dispersiveness are also discussed. 

Among the stable systems, Poisson and Liapunov stability 
will be discussed as well as non-wandering systems. Poisson 
instability, Lagrange instability, and complete instability 
will also be investigated. Topics in dispersiveness include 
parallelizable and dispersive systems. The relationship 
between several of these systems will be shown. One such 
relationship to be shown is that a parallelizable system is 
dispersive, completely unstable, Lagrange unstable and 
Poisson unstable. 

Several important and revealing theorems will be pre- 
sented. Among these are the theorems involving the Liapunov 
function which allows analysis of a system as to stability 
or instability without actually solving the system. 



The book Stability Theory of Dynamical Systems by 
Bhatia and Szego will be used extensively throughout this 
thesis and its bibliography contains an excellent list of 
references . 



i 
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I. TOPOLOGICAL PROPERTIES 



The purpose of this section is to expose the reader to 
all of the topological properties needed to understand the 
material presented in this thesis. It is important that 
the reader develop a firm understanding and familiarity 
with these concepts to facilitate in understanding of the 
material to follow in later sections. 

We begin with the concept of distance. 

If X is a non-empty set, a metric or distance function, 
on X is a real-valued function d of ordered pairs of elements 
of X which satisfy three conditions. First, d(x,y) j> 0, 
and d(x,y) = 0 if and only if x = y. Second, d(x,y) = d(y,x). 
This is called the symmetric property . Finally, 
d(x,y) £ d(x,z) + d(z,y). This property is called the 
triangle inequality . A metric space X is then a set X with 
a- metric defined on it. 

If Xq is any point of a metric space and r is a positive 
real number, the open sphere S(xQ,r) with center Xq and 
radius r is the subset of X defined by S(xQ,r) = {x:d(x,XQ)<r} . 
Now a subset G of X is said to be open if and only if it is 
the union of open spheres. Since Xq is clearly an element 
of S(xQ,r), we say that S(xQ,r) is a (spherical) 
neighborhood of x^. 

If A is a subset of a metric space X, a point x in X is 
called an accumulation point of A if each open sphere 
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centered on x contains at least one point of A different 
from x. The essential idea here is that the points of A 
different from x get arbitrarily close to x. Then a subset 
F of a metric space X is called a closed set if it contains 
all its accumulation points. Another characterization of 
a closed set is that a set F C X is. closed if and only if 
its complement in X is an open set . Then by definition, 
the whole space X and the empty set 0 are open sets. Then 
since the empty set is the complement of the whole space, 
we have by our previous observation that the whole space X 
and the empty set 0 are also closed sets. 

Now, in any metric space, the union of any collection 
of open sets is itself an open set and the finite inter- 
section of open sets is an open set. Taking set complements 
it follows that any intersection of closed sets is a closed 
set and the finite union of closed sets is a closed set. 

If Xq is a point in the metric space X and'r is a non- 
negative real number, then the closed sphere S[xq,t] with 
center Xq and radius r is the subset of X defined by 
S[x Q ,r] = {x:d(x,x Q ) < r] . 

Let X be an arbitrary metric space and let A be a subset 
of X. A point in A is called an interior point of A if it 
is the center of some open sphere contained in A; and the 
interior of A, denoted by int^A, is the set of all its 
interior points. Symbolically, we have 

int^A = (x:x e A and S(x,r) C A for some r} . Since the 
interior of a set is a union of open spheres, it follows 
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that int^A is an open set. Hence we have the characterization 
that a set A is_ open if and only if A = int^A. We often 
write int A in place of int^A whenever no confusion will 
result . 

If X is a metric space and A is a subset of X, the 
closure of A in X, denoted by cl^A, is the union of A and 
the set of all its accumulation points. Intuitively, cl^A 
is A itself together with all other points in X arbitrarily 
close to A. Since it is clear that cl^A contains all its 
accumulation points, cl^A is a closed set. Thus we have 
an additional characterization of a closed set : AC X is 

closed if and only if A = cl^A. We often write cl A in 
place of cl^A. 

Again let X be a metric space and let A be a subset 
of X. We will use the notation X-A to denote the complement 
of A in X. A point in X is called a boundary point of A 
if each open sphere centered on this point intersects both 
A and its complement X-A. Now the boundary of A, denoted 
bdyA, is the set of all the boundary points of A. It is 
easily shown that bdy A = cl^A 0 cl^(X-A) . Moreover, the 
boundary of A is a closed set, and a set AC X is closed 
if and only if it contains all of its boundary points . 

Every metric space has, in addition to the properties 

already stated, two important "separation" properties. 

First the Hausdorf f property: if x and y are two distinct 

points in X, then there exist disjoint open sets U and U 

a y 
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such that x £ U„ and y e U„. Second, the normal space 

x y 

property: if A and B are disjoint closed sets of X, then 

there exist disjoint open sets U A and Ug such that A C U A 
and B d Ug. 

We now use these basic definitions to describe more 
advanced ideas in topology. We are interested in character- 
izing "compactness", "completeness" and "connectedness". 

Let X be any metric space. A class (G^) of open subsets 

of X is said to be an open cover of X if each point in X 

belongs to at least one G. , that is U G. = X. A subclass 

1 i 1 

of an open cover which is itself an open cover is called a 
subcover . Then a compact space is defined to be a metric 
space in which every open cover has a finite subcover. This 
becomes a familiar definition when one recalls the Heine - 
Borel Theorem which states that if [a,b] is_ a closed interval 
of real numbers and if (G, } , x e [a,b] is an open cover of 
[a,b], then there exists a finite number of the G„'s whose 
union contains [a,b]. This is just the statement that 
closed intervals on the real line are compact sets. 

A much used property of a compact space is that any 
closed subspace of a compact space is itself compact . Another 
useful property is that any infinite sequence in a compact 
set contains a subsequence which converges to some point 
in the compact set. This also becomes a familiar property 
when one recalls the Bolzano - Weierstrass Theorem : Every 

bounded infinite set A of real numbers has at least one point 
of accumulation . This theorem is just an application of 
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this property applied to the real line. These last two 
properties are very important and will be used extensively 
throughout this thesis. 

To define the notion of a complete metric space we need 
the concept of convergence. If {x n > is a sequence in X, 
we say that {x n > is convergent to the point x in X if given 
any e > 0, there exists a positive integer N such that 
d(x n> x) < e whenever n > N. Then in particular we have 
d(x n ,x m ) < e whenever n,m > N. A sequence with this latter 
property is called a Cauchy sequence . Now, not every Cauchy 
sequence is a convergent sequence. For example, consider 
the following sequence in the metric space of rational 
numbers, (3., 3.1, 3.1*1, 3.1*11, 3.1*115, ...}. Clearly, this 
is a Cauchy sequence but it is not convergent because it 
converges to ir; but v is not a rational number. Motivated 
by this example we say that a metric space is complete if 
every Cauchy sequence in X is convergent. Completeness is 
related to compactness in that a set is compact if and only 
if it is complete - and totally bounded. By totally bounded 
we mean that the space X is the union of a finite number of 
open spheres of radius less than r for each positive r. 

Finally, we define the concept of connectedness. A 
connected space is a metric space X which cannot be repre- 
sented as the union of two disjoint non-empty open sets. 

If X = A U B where A and B are disjoint non-empty open sets, 
then X-A = B and X-B = A so A and B are also closed sets. 
Thus a space X is_ connected if and only if the only sets 
which are both open and closed are the whole space X and the 
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empty set 0. This property is very important and it can be 
shown that the Intermediate value theorem in calculus is 
equivalent to this property for intervals on the real line. 

We now turn our attention to the concept of a continuous 
function . This concept has several different characteriza- 
tions and we begin with the basic definition that looks like 
the definition used in elementary calculus. 

Let f be a function from a metric space X with metric d 
to a metric space Y with metric d'. Then f is said to be 
continuous at Xq e X, if for each e > 0, there exists a 
6 > 0 such that for x e X, we have d ' ( f ( x) , f ( X q ) ) < e 
whenever dCxjX^) < 6. The function f is said to be continuous 
if it is continuous at each point of X. The definition we 
saw in calculus was for real-valued functions so the metrics 
just involved the usual absolute value for the real line. 

Several other characterizations of a continuous function 
will be needed in this thesis and are stated here. A function 
is continuous if and only if the inverse image of each open 
( closed ) set is an open ( closed ) set . Also a function f is 
continuous if and only if for each sequence (x n ) which 
converges to x , the sequence { f ( x ) } converges to f ( x ) . 

We now state some important properties of continuous 
functions . 

A continuous real - valued function f 
defined on a compact set A attains its 
maximum and minimum values on A in the 
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following sense : if a = inf{f(x):x e A} 

and b = sup{f(x):x e A}, then there exists 
x^,x^ e A such that f(x^) = a, and 
f(x 2 ) = b. 

The continuous image of a compact ( connected ) 
set is compact ( connected ) ♦ 



Next, we define the concept of a homeomorphism. A 
homeomorphism is a one-to-one continuous function from one 
metric space onto another metric space such that the inverse 
of this function is also continuous. (A one -to- one function 
is a function which assigns different images to different 
elements in its domain. A function f from a metric space 
X to a metric space Y is onto if the image of X under f is 
the entire space Y. If a function is one-to-one and onto, 
then its inverse is well defined.) 

We have already stated the definition of a compact space. 
A closely related idea is that of a locally compact space. 

As its name suggests, a metric space X is locally compact 
if and only if each point of X has at least one compact 
neighborhood. Hence, a compact space is automatically 
locally compact. 

One of the most important properties of a locally compact 
space is that through the addition of one point, usually 
referred to as the ideal point, the space can be extended 
to a compact space. This extension is known as the one-point 
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compactification . The one - point compact ificat ion of a 
metric space X is the space X' = X U {to} where the open sets 
in X' are all the open sets in X together with all subsets 
U of X* * such that X’-U is a compact subset of X. In other 
words, a set U is open in X' if and only if U f\ X is open 
and whenever to e U, X-U is compact. As this is an important 
fact, we give a formal proof. We must show that X* is compact 
and that unions and finite intersections of open sets are 
open. 

THEOREM : The one - point compactification X' of a metric 

space X is_ compact . 

PROOF: Finite intersections and arbitrary unions of open 

sets in X* intersect X in open sets. If w is a member of 
the intersection of two open subsets of X', then the comple- 

* ment of the intersection is the union of two closed compact 
subsets of X and is therefore closed and compact. Hence 
the intersection of two open subsets is open in X'. If to 
belongs to the union of the members of a family of open 
subsets of X', then to belongs to some member U of the family. 
Moreover the complement of the union is a closed subset of 
the compact set X-U and is therefore closed and compact. 

Hence the union of open sets is open in X'. Consequently, 
arbitrary unions and finite intersections of open sets in 
X* are open. Next, let be an open covering of X*. Then 
w is a member of some U in U and X-U is compact. Hence 
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there is a finite subcover of W- which covers X-U and thus 
this finite subcover together with U is a finite subcover 
of X'. Thus X’ is compact. // 

Our final topic in this section of topological prereq- 
uisites is the concept of product spaces. A finite product 
space is the cartesian product of two metric spaces. In 
this thesis, one of the metric spaces will usually be the 
real line R; consequently if X is a metric space, we have 
the product space to be XxR. Hence any point in the product 
space will be written (x,t) where x e X and t c R. Now the 
open sets in the product space are just the cartesian products 
of open sets in the metric spaces. That is if U is an open 
set in X and V is an open set in R, then UxV is an open set 
in XxR, and every open set in XxR is of that form. If d is 
the metric for X, we define the metric d* on XxR in the 
following manner. If x,y e XxR, then there are elements 
a,b e X and r,s e R such that x = (a,r) and y = (b,s). Then 
d’(x,y) = d(a,b) + |r - s|. The spaces X and R are called 
the coordinate spaces of the product and the functions 
and Pp, which carry the point (x,t) of XxR into x and t 
respectively, are called the projections into the coordinate 
spaces. It is easily shown that the projection mappings 
are continuous functions and that their inverses are also 
continuous functions. Moreover f:Y X is continuous if 
and only if each of the compositions P^of and Py°f is 
continuous, for any metric space Y. 






A function f = f(x,t) defined on the product space is 
continuous in x if and only if for each t the function f( ,t) 
whose value at x is f(x,t), is continuous. Similarly, 
f(x,t) is continuous in t if and only if for each x e X, 
the function f(x, ) such that f(x, )(t) = f(x,t), is con- 
tinuous. Finally if f is continuous on the product space 
then f is continuous in both x and t by our previous remarks. 
However, the converse is not, in general, true. 
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II. DYNAMICAL SYSTEMS 



We seek to motivate the concept of dynamical systems 
by some examples of differential equations. But first we 
need the existence and uniqueness theorem associated with 
such systems. 

We consider systems of first-order differential equations 
written in the form 

dx, 

dt~ = f' 1 (t;x 1 , . . . ,x n ) 



dx 

-nr^- = f (t :x, , . . . ,x ) . 
dt n * l 5 5 n 

Such systems frequently arise out of physical problems where 
the independent variable t represents time. Then a solution 
of the above system describes the "trajectory" of a particle 
moving in n-space. To see this more clearly 3 we look at 
the following two examples. 

•3 

Consider the following system defined in R : 




-sint 



dx 



2 



dt 



cost 



1 . 



dx 

dF 



3 _ 



Then the general solution to this system is given by 



x^ = cost + 



x 2 = sint + c 2 



Xj = t + Cj 



where the c's are all constants of integration. 

Suppose we specify that this solution passes through the 
point x^ = 1, x 2 = 0, and x^ = 0 at time t = 0. Then the 
solution is given by 



x 1 = cost 



x 2 = sint 



x 3 = t. 



Such a solution describes the "trajectory" sketched below. 
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It is clear that at any time, we may find the exact 
position of the particle. The same is true of the following 
example . 

Consider the system defined in R given by 




dx 2 

dF" 



2t 



dx? 

dtT 



t 



2 



where we specify that its solution passes through the point 
X 1 = x 2 = and x 3 = 0 at time t = 0, Then this system 
has the solution given by 



x 



1 



2t 
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Its "trajectory" is sketched below. 

*3 




Thus an initial value problem for such a system consists 
of selecting an initial time t = tg and a point Xq in R n , 
and then seeking a solution whose graph passes -through Xq 
at time t = tQ. 

By setting X(t) = (x^t ) , . . . ,x n (t ) ) , 

f(t,x) = (f 1 (t,X),. . . ,f n (t,X)), 

this system may be rewritten as 

gjr = f(t,X), where = (x^^ ' (t ) , . . . ,x fi ' (t ) ) . 

We next want to examine conditions under which unique 
solutions exist for an initial value problem of the system. 
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Suppose that f = f(t,X) is bounded and continuous in a region 

n+ n 

D of R n . Furthermore, assume that f satisfies a Lipschit z 
condition in D; that is there exists a non-negative real 
number b such that whenever (t,X^) and (t,X2) are any two 
points in R, d(f(t,X^), f(t,X2)) £bd(X^,X2) where d is 
the ordinary Euclidean metric. Then we have the following 
theorem. 



Let f = f(t,X) be continuous , bounded, and 
satisfy a Lipschitz condition in a region D 
of R n+ \ and let (t^^^) be any point in D. 
Then the initial value problem 



dX 

dt 



f(t,X) 



X(t Q ) = X 



0 



has a unique solution on ari interval I of R 
and the graph of the solution lies in D. 



We will be interested mainly in systems defined on all 
of R. Then by ensuring that f is continuous and bounded on 
all of R n and satisfies a global Lipschitz condition, we 
will have existence and uniqueness of solutions to the above 
system. 

One elegant and simple proof of the above theorem requires 
a knowledge of contraction mappings and some theory of 
fixed points and will not be presented here [Kreider, Kuller, 
Ostberg, p. 385]. 



20 



A classical proof may be found in Theory of Ordinary 
Differential Equations by Coddington and Levinson. 

Of critical importance is the question: what happens 

to the solution if small changes are made in the function f 
or in the initial conditions t Q and X Q ? In most physical 
problems whose solutions are obtained from such a system, 
experimental error and estimations of physical constants 
may affect the initial conditions, and the presence of 
variable parameters, such as temperature or density, may 
affect f. We would hope that such small variations would 
not cause a drastic variation in the solution. It can be 
shown that these unique solutions guaranteed by the theorem 
depend continuously on f and on t^ and Xq for a certain 
amount of time; that is small changes in f or in t Q and Xq 
produce small changes in the solution. More precisely. 

Let f = f(t,X) be_ continuous , bounded , and 
satisfy a Lipschit z condition in a domain D 
of the (n + 1) - dimensional (t,X) space , and 
suppose X is_ a solution of the system on an 
interval I : a <_ t £ b . There exist s a 6 > 0 
such that for any (^,Xq) where a < t^ < b, 
and dCXQjXCt^ ) ) < <5 there exists a unique 
solution Y on I with Y(t 2 3 t^,X 0 ) = X Q where 
a < t 2 < b. Moreover , Y is_ continuous in D. 

(Proofs of the continuity of solutions may be found in 
Kreider, Kuller, Ostberg, p. 390 and Coddington, Levinson, 
p. 22.) 
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What happens to these solutions over an extended period 
of time involves the theory of stability which is the primary 
focus of this thesis. 

We will be interested mainly in first-order autonomous 
systems , that is, "time independent" systems. Such a system 
may be written in the form 



X = f(X) where X denotes ^ 
x(t 0 ) = X 0 . 

The following examples provide some insight into such 
autonomous systems. 



EXAMPLE 2.1 Consider a simple pendulum of mass m swinging 
without friction on a weightless rod of length 1. By 
Newton’s second law, the equation of motion of the pendulum 
is given by 



-mgsin 6 = ml 




0 = -g/1 sin 0 
2 

By setting g/1 = k and 0 = co we 
obtain the following system: 
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0 = to 



* 2 

to = -k sin 0. 

It can be shown [Kreider , Kuller, Ostberg, p. 397] that 
the general solution to this system is 

to 2 = 2k 2 cos 0 + c 

2 

where c -2k is a constant depending on the initial 
conditions . 



Graphing several of these solutions in the 0to-plane we 
have the following diagram. 

w 




To see how this diagram was obtained, we see that to is 

2 

the angular velocity. Suppose c = -2k . Then when 0=0, 
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a) = 0 and hence 0 and w are both zero. Thus the system is 

not changing with time. The point (^) is the complete 

trajectory for this case. Physically this shows that if 

the mass is started at the equilibrium position with no 

angular velocity, the mass will remain at the equilibrium 

2 

position indefinitely. Suppose c = 2k and the system is 

» • 

started at 0 = tt. Then again w = 0 and hence 0 and w are 
both unchanging with time and the point (q) describes the 
complete trajectory. Physically this means that if the 
mass is balanced vertically above the equilibrium point 
with no angular velocity, it will also remain in that position 
indefinitely. Clearly the same analysis can be made for 
6 = ±mr which is exactly what we would expect. If this 
system is started at any other angle 0, it will have an 
initial velocity which will cause the mass to approach the 
angle 0 =tt or 0 = -it and it will get arbitrarily close 
but will never actually reach that point. Neither will this 
mass swing back away from this point but rather it will 
continue to approach the angle 0 = tt or 0 = -tt as time goes 
to infinity. This is the physical meaning of the trajectories 

between the points ) and (q). 

2 2 

Now suppose -2k < c < 2k . From the previous remarks 

it is clear that the mass will have a velocity at 0 = 0 but 
its velocity is not great enough to permit it to approach 
the angle 0 = ir or 0 = -ir. We would expect in this case 
that the mass would swing back and forth across the equili- 
brium point always reaching the same angle on either side 
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of the equilibrium point. This is the physical interpre- 
tation of the elliptical trajectories shown in the diagram. 

2 2 

Finally, suppose c > 2k . Then w is always positive 
which means that w is always positive or always negative. 

If to is always positive, then 0 is always positive and 
hence 0 is always increasing. A similar statement follows 
if w is always negative. Physically this is the case when 
the mass has been given enough initial velocity to swing 
completely over the top position. In a frictionless envir- 
onment we would expect this behavior to continue indefinitely. 
Such is the physical interpretation of the top and bottom 
trajectories in the diagram. 

To see the direction of a particle along any trajectory, 
pick an angle 0. For 0 < 0 < tt, we see that w < 0 so to is 
decreasing. For -tt < 0 < 0, we see that to > 0 so to is 
increasing. This completes our analysis of the diagram. 

It is noted that in this example the solution is not 
given as a function of time. To obtain such would mean the 
difficult task of solving a non-linear differential equation. 
This is not any more informative since the movement of a 
point along any trajectory with increasing time can easily 
be determined from the system describing the motion as we 
have just observed. 

A more elementary approach to the problem is by observing 
that for small angles 0, sin 0 is approximated by 0. Thus 
our system would be given by 
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9 = 0 ) 



* 2 

o) = -k 0. 

It can be shown that the solution to this system is 

9 = Asin kt + Bcos kt 
o) = Akcos kt - Bksin kt 

where A and B are constants depending on the initial conditions 

[Spiegel, p. 210]. It is seen that this solution describes 

simple harmonic motion and closely approximates our trajec- 

2 2 

tories for small c where -2k < c < 2k . 

EXAMPLE 2.2 We consider the following vibrating spring 
problem. 




Wy 

g 



+ by + ky = 0 
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where W is the weight of the mass, g is the force of gravity, 
b is the damping constant and k is the proportionality 
constant of the spring. Suppose W = 6 lbs, g = 32 lbs, 
k = 12 and b > 0 is left as an undetermined constant for 
now. 

Then it can be shown [Spiegel, p. 195] that the general 
solution for y is 



where the constants A and B depend on the initial conditions. 

If we specify the initial conditions to be y(0) = 1/3 
and y(0) = 0, we have the solution, 



To make this example more explicit we will consider two 
cases for b. First an undamped system with b = 0 and then 
a damped system with b = 1.5* With b = 0, we have 



y - e 



^ (Acos ( - j J 9 - b^ t) + Bsin ( ^- J 9 - b^ t)). 




y = ^ cos 8t 



With b = 1.5, we have 



y = ^- e“^ t (3cos 4 t + ^ sin 4 J"3 t ) 



By use of a trigonometric identity, this last equation may 



be written 



2 J~3 



y = — j - - e ^ sin (4 J~3 t + 



Now that we have the solutions, we change this second- 
order problem into a system of first-order equations. Let 
X = ( X ) and y = x. Then for b = 0 

y 

iM y 

. w J 

X = f(X) = ( ) 

x 



x(°) = q /3 ) 



For notation, let the initial point (^ 73 ) = X Q . Then the 
solution to this system is 



X(X Q ,t) 



-i-sin 8t 

( 3 ) 

^cos 8 t 



For b = 1.5, we have 



X = f(X) = 



-bgx - kgy 
W “ W 



x 



x(o) = q /3 ) 



and the solution to this system is 
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x(x Q ,t) = ' 3 



^cos - |-J3e ^sinC*! Jlt+^O 



|-i3e itt sin( i l43t + j) 



It is clear that as t t «>, X(Xg,t) -*• (g) which we would 
intuitively expect from a damped spring. 

We sketch these solutions on the x,y-plane. 





FIGURE 2.2 



For the undamped problem we see that its trajectory 
describes simple harmonic motion. It continues to oscillate 
indefinitely with the same maximum displacement from the 
equilibrium position. On the other hand, the displacement 
and its rate of change both decrease with increasing time 
for the damped problem. 

Since we have two different systems, it is necessary to 
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plot their solutions on two different planes. This helps 
to avoid any misunderstanding in our future work. 

EXAMPLE 2.3 Consider the second-order differential equation 
y + y = 0, where y is considered as a function of time. 

By setting x = y, this equation may be converted to the 
system 



x = -y 



Now considering X = ( ) with the initial condition that at 

y 

t = 0, the solution passes through the point (q), we have 
the following initial value problem 



X = f(X) 



X(0) = (J). 



where f is the linear function defined by 




(”^). A simple check shows that 

X 



this system has the solution 




If we let X Q denote the point (q), then Xq denotes the 
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initial condition. Further since the solution depends on 
this initial condition, the solution X(t) is more properly 
stated X(X Q ,t) = ( I ? ® £). 



There are three observations we wish to make about this 
solution. First, we see that X(X Q ,0) = X Q . Next, 
X(X(X 0 ,t 2 ),t 1 ) = X(X 0 ,t 1 + t 2 ). 

cos t p 

This is true since X(X(X Q , t 2 ) ,t ) = X( ( ± t ),t 1 ). 



The term on the right is the solution of the system which 

cos t 2 

has the initial condition Yq = ( s ^ n ^ )• Its solution is 



X(Y 0 ,t) 

X(Y 0 ,t 1 ) 



cos t 0 cos t - sin t 0 sin t 

( ^ ) 
cos t 2 sin t + sin t 2 cos t 

cos t 2 cos t^ - sin t 2 sin t-^ 
^cos t 2 sin t^ + sin t 2 cos t^ 



Hence 



) . Using a 



trigonometric identity we have 



X(Y Q3 t 1 ) = ( 



cos 

sin 



(t l + 

(t x + t 2 r 



Thus 



X(X(X 0 ,t 2 ),t 1 J = X(Y 0 ,t 1 ) = X(X 0 ,t 1 + t 2 ). Finally X(X Q ,t) 
is continuous because both of its components are continuous. 

It is clear that the graph of the solution is a circle 
in the x,y-plane, and if we consider the graph to be the 
trajectory of the point Xq , that point moves counter-clockwise 
with increasing time. This can be seen by choosing some 
positive value for x. Then y is positive so y is increasing. 
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'I 




FIGURE 2.3 

£ 

If we alter the initial condition to X^ = ( fe ) , the 

a cos t - b sin t 

solution is X(X Q ,t) = ( sin t + b cos t ) and its trajectory 

is a circle passing through the point X^. 

Now referring to Example 2.3 we see that the solution 

2 

X(X Q ,t) is a map from the product space R xR into the space 
2 2 

R . We note also that R is a metric space. The three 

observations made about the solution are precisely the 

properties defining a dynamical system. The fact that the 

system in Example 2. 3 does define a dynamical system will 

be made clear by the following definition where in this case, 

the map tt is defined by ir(x,t) = X(x,t) where x denotes a 
2 

point in R . In all further work, X denotes an arbitrary 
metric space. 

DEFINITION 2.^ A dynamical system on X is the triplet 
(X,R,ir), where tt is a map from the product space X*R into 
the space X satisfying the following axioms: 

1. For every x in X, tt(x, 0) = x. 
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2. For every x in X and t^,t 2 in R, 
iTC-irCxjt^) ,t 2 ) = TrCxjt^ + t 2 ), and 

3. The map tt is continuous. 

The space X and the map tt are called the phase space and 
the phase map , respectively. 

From this point on, whenever no confusion will result, 
the symbol tt will be suppressed. Thus the image under tt 
of the point (x,t) will be written more simply as xt . The 
first two axioms then become 

1. For every x in X, xO = x, 

2. For every x in X and t^,tg 
xt ^ ( 1 2 ) = x ( t ^ + t 2 ) . 

Similarly, if M <Z X and A R, MA is the set 
{xt : x e M and t e A} . 

The phase map determines two other maps when either x or 
t is fixed. For fixed t in R the map tt :X x defined by 
tt^(x) = xt is called a transition ; for fixed x in X the 
map tt..:R -* X defined by Tr..(t) = xt is called a motion 

X X ' 

through x. 

THEOREM 2.5 For each t iri R, tt^ is_ a homeomorphism on X onto 
itself . 

PROOF For any t in R, Tr t is clearly continuous because tt 
is continuous. tt^ is injective since xt = yt implies x = y 
because x = xO = x(t - t) = xt(-t) = yt(-t) = y(t - t) = yO 
=y. To see that tt^ is surjective, let y e X be arbitrary. 
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Then ii t (x) = y where x = y(-t). Finally to show that tt^ 
has a continuous inverse, it is sufficient to show that 

is the inverse of tt^. To see this we note that if t 

s t s 

and it are two transitions, the composition tt o-rr is the 

t s 

transition tt : for any x in X, 



7 T t o Tr S (x) = TT t (lT S (x)) = TT^Cxs) = Xs(t) = x(s+t) = 1T S + t (x). 

It is also clear that n ® is the identity transition because 

for any x in X, tt^(x) = xO = x. Now since 

“t t t t 0 —t/ 

7 T ° it = 7 T = if , the transition ir is the inverse of 

t 

IT . 

REMARK 2.6 For any x in X and [a,b] C R, the set x[a,b] 
is compact and connected. 

PROOF The set {x} is compact and connected in X and [a,b] 
is compact and connected in R. Thus {x}x[a,b] -is compact 
and connected. Then since x[a,b] is the continuous image 
of this set, x[a,b] is compact and connected. // 

REMARK 2.7 Observe that the transitions tt^ , t in R, form a 
commutative group with the group operation being the compo- 
sition of transitions. 

PROOF The closure property is trivial for if t,s are in R, 

+- s t s 

t + s is in R and tt ott = tt as we have already observed. 
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Associative : 



t , s v, 

TT o (. IT OTT ) 



t + S V 
= IT °TT 



t s+v 

TT ott 



/ t Sv V 
(iT ott Jott 



t + S + V 
TT 



Commutative: TT t o 7 T s (x) = TT t (ir S (x)) = tt^xs) = (xs)t 

= x(s+t) = x(t+s) = (xt)s = ir S (xt) 

= TT S (tT^(x)) = TT S OTT^(x) 

Identity: The identity is tt^ since tt^o-tt^ = tt^ + ^ 

t 0+t 0 t 

= IT = TT = TT °TT 



t — 

Inverses: The inverse of tt is tt since 

t -t 0 

TT oTT - TT . 



We now give a precise definition of what we have been 
calling trajectories. 

2.7 DEFINITION The maps y(x), y + (x), and y~(x) from X into 

X 

2 are defined, for any x in X, 

y(x) = {xt: t e R}, 

Y + (x) ='{xt: t e R + ), 

Y~(x) = {xt: t e R - }. 
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For any x in X, the sets y(x), y + (x), and y”(x) are called 
respectively the traj ectory , the positive semi-traj ectory 
and the negative semi - traj ectory . Note that for any x in X, 
y(x) = xR. 

If two trajectories and Y 2 have a point in common, 
then the trajectories are identical. This fact is guaranteed 
by the existence and uniqueness theorem. For if x e y-j_ 
and x e Y 2 > this theorem states that there exists a unique 
solution passing through the point x. Hence y-j_ = ^ 2 * This 
means that two different trajectories can never cross each 
other. This is an important fact and is very useful in 
diagramming dynamical systems. 

Another example of a dynamical system which will prove 
to be useful is a dynamical system defined on a torus. 
Consider a differential system in the plane 

X = f(X) 



where f satisfies the hypotheses of the existence and unique- 

f. 



ness theorem. Set X = ( x ) and f = ( f ' L ) and assume that the 

y x 2 

functions f^ and are periodic with period 1 in each of 
the variables x and y. Then 



f^Xjy) = f i (x+l,y) = f^ (x ,y+l) , i = 1,2 

V/e take the square {(x,y): 0. £ x < 1, 0 £ y < 1} and 
identify the sides x = 0 and x = 1 and the sides y = 0 and 
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y = 1. Then the points (0,0), (0,1), (1,0) and (1,1) are 
identified. Then we project the plane onto the torus by 
identifying any point (x,y) in the plane with a point (x,y) 
of the torus where x = x (mod 1) and y = y (mod 1). Then 
the trajectories in the plane projected onto the torus yield 
trajectories of a dynamical system on the torus. 




It is clear from Example 2.3 that if we pick any point 
x on the trajectory, xt is also on the trajectory for any t 
in R. In Example 2.2 if we pick a point x on any trajectory 
and consider the positive semi-trajectory, then xt belongs 
to y + (x) for any t e R + . The same results can be seen for 
the negative semi-trajectory. This observation leads us 
naturally into the following concept of invariance. 

2.8 DEFINITION A set M C X is called invariant whenever 

xt e M for all x e M and t e R. 

It is called positively invariant whenever this holds with 
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R replaced by R + and Is called negatively Invariant if it 
holds with R replaced by R~. 

From this definition we see that a set H C X is invari- 
ant if and only if it is both positively and negatively 
invariant for if M is invariant, we let x be an arbitrary 

element of M. Then xt e M for all t e R and in particular 

* 4 * — 

for t e R , xt e M and for t e R , xt e M. Thus M is both 

positively and negatively invariant. If M is both positively 

+ 

and negatively invariant, then for any t e R , xt e M and 
for any t e R , xt e M. Thus xt e M for all t e R and M 
is invariant. 

This concept leads us to some interesting and important 
properties associated with such sets. 

THEOREM 2.9 Let {M^} be a collection of positively invari- 
ant , negatively invariant or invariant subsets of X. Then 
their intersection and their union have the same property . 

A set M C X is_ positively invariant if and only if the set 
X-M i_s negatively invariant . M is invariant if and only 
if X-M is_ invariant . 

PROOF Let the sets M^ be positively invariant. Let A =UM^ 
and B = Am^. For any x e A, x e for some i. Thus 

xt e M^ for all t e R + since M^ is positively invariant. 

Then xt e A for all t e R + . Hence A is positively invariant. 
Let x e B. Then x e M^ for all i. Then since each M^ is 
positively invariant, xt e M^ for all t e R + and all i. 
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Thus xt e B for all t e R + and B Is positively Invariant. 

The proofs for the negatively invariant and invariant cases 
are entirely analogous. Let M be positively invariant. 

Let x e X-M and t e R*" and suppose xt e M. Since -t e R + , 
xt(-t) e M. But xt(-t) = xO = x which implies x e M and 
we have a contradiction. Thus xt e X-M and X-M is negatively 
invariant. The proofs of the converse and the third part 
are entirely similar. // 

THEOREM 2.10 Let M C X be_ positively invariant , negatively 
invariant or invariant . Then the following statements are 
true . 

1. The closure of M, cl M, and its interior , int M, 
have the same property , 

2. Each of its components has the same property . 

If M is_ invariant , then so is its boundary bdy M. 

PROOF 1. Consider the case of invariance. Let x e cl M 

and t e R. Then there is a sequence (x^) in M such that 

x •*- x. By the invariance of M, x t e M for all t e R 
n J n 

and each n. Then by continuity we have that x t ■+ xt. Thus 
xt e cl M and cl M is positively invariant. The proofs for 
the negatively invariant and positively invariant cases 
are entirely analogous. 

Assume M is positively invariant. Then X-M and 
hence cl (X-M) are negatively invariant. Consequently 
int M = X - cl (X-M) is positively invariant. The proof for 
M negatively invariant is entirely similar. 



39 



2. Let M be invariant. Now M = UM^ where the M^'s 
are disjoint components of M. Suppose M i is not invariant 
for some i. Then there exists x e M.^ and t e R such that 
xt £ M^. Nov; IYL is a maximal connected subset of M and 
x[0,t] is a connected set. Thus if xt e M^ for some j, 

M. U M. is a connected set and M. and M. are not disjoint. 

1 J 1 J 

Thus xt M. for all j v;hich implies that xt i M. This 

J 

contradicts the invariance of M and hence each component of 
M is invariant. For the converse let each component of M 
be invariant. Then since M = and each component is 

invariant, it follows that M is invariant. The proofs of 
the other cases are similar. 

Suppose M is invariant, then so is X-M. Consequently 
cl M and cl (X-M) are invariant. Hence bdy M = cl M A cl (X-M) 
is invariant and int M = X - cl (X-M) is invariant. // 

For the following theorem we need to define y(M), y + (M), 
and y (M) where M C. X. We define these as we would expect; 
that is y (M) = {Mt : t e R} , y + (M) = {Mt: t e R + } and 
similarly for y~(M). 

THEOREM 2.11 For any x e X, the sets y(x), y + (x), and 
y (x) are , respectively , invariant , positively invariant 
and negatively invariant . A set M C. X is invariant , posi - 
tively invariant , or negatively invariant if and only if , 
respectively , y(M) = M, y + (M) = M, or y (M) = M. 
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PROOF Let x^ e y”(x) for e R _ and let t 2 e R”. Then 
xt^(t 2 ) = x(t^ + t 2 ) = xt^ for t^ = t^ + t 2< Since 
t 1 ,t 2 e R , we have t^ e R*~ and xt^ e y*"(x). Hence y“(x) 
is negatively invariant. The proof of the other cases is 
entirely analogous. For the second part, let M be negatively 
invariant. Nov/ y**(M) = (Mt : t e R - }. In particular for 
t = 0, Mt = M and MC y”(M). Suppose y”(M) <£ M. Then there 
exists a t^ e R such that Mt^ M. Then there exists 
x e M such that xt^ 0 M which contradicts the negative 
invariance of M. Thus y**(M) C M and y”(M) = M. Conversely, 
suppose y (M) = M. Then (Mt : t e R~} = M and Mt C m for 
all t e R _ . Thus M is negatively invariant. The proofs of 
the other cases are similar. 

The following remarks are useful in understanding this 
material and prove some useful results. 

REMARKS 2.12 

1. The sets X and 0 are invariant. 

PROOF Let x e X. Then xt e X for all t e R since xt 0 0. 
Thus X is invariant. Since 0 = X-X, it follows that 0 is 
invariant . 



2. A set M C X is invariant, positively invariant, 
or negatively invariant if and only if for each x e M, 
respectively, y(x) C M, y + (x) CM, or y (x) C M. 
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PROOF Let M be invariant. Then by Theorem 2.11 y(M) = M. 
Thus {Mt: t e R} = M. Hence for any x e M, xt e M for all 
t e R. But {xt: t e R} = y(x). Thus y(x) C M. Conversely, 
suppose y(x) c. M for each x e M. Then {xt: t e R} c M and 
xt e M for all t e R and all x e M. Hence M is invariant. 

The proofs of the other cases are similar. 

nT r 

Consider once again Example 2.1. For the points ) 

' ' 2 

n = ±1, ±2, ... we see that 0 = w = 0 and to = -k sin 0 = 0. 

I1TT 

This shows that each x = (q ) does not change position with 
increasing or decreasing time. Such points are called 
critical or rest points. We now state the formal definition. 

DEFINITION 2.13 A point x e X is said to be a critical 
point if x = xt for all t e R. 

REMARK 2.14 If x is critical and y f x, then x does not 
belong to the trajectory y(y). 

PROOF Suppose x is critical and y ¥ x. Then {x} = y(x) 
and since y(x) is invariant, {x} is invariant. Thus 
X-{x) is invariant. Since y j* x, y e X-{x) so y(y) C X-{x). 
Thus x t y(y). // 

With this remark. Figure 2.1 should more properly be 



shown as follows : 



The next theorem contains several characterizations of 



critical points but first the following lemma will be found 
useful . 

LEMMA 2.15 If x e X and x = xt for some t e R, then 
x = x(nt) for all ' integers n. 

PROOF If x = xt j then x(-t) = xt(-t) = x(t - t) - x so we 
need only consider positive integers. We complete the 
proof by induction. Clearly if x = xt for some t e R + , 
then x = x(lt). Assume x = x(nt). Then 

x=xt(nt)=x(t+nt)=x((n+l)t). // 

THEOREM 2.16 Let x e X. Then the following are equivalent . 

1. The point x is critical 

2. The singleton (x) is_ y(x) 
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3. 


The singleton 


(x) 


is 


Y + (x) 


4. 


The singleton 


{ x) 


is 


Y"(x) 


5. 


The singleton 


(x> 


is 


x[a,b] for some a < b 


6. 


There is a sequence 


(t }, t > 0, t ->-0 with 
n 5 n * n 



x = xt for each n. 



PROOF The equivalence of statements 1. with 2., 3. > and 4. 

are trivial. We next prove the equivalence of statemnts 1. 

and 6. If x is critical then statement 6 follows trivially 

since x = xt for all t e R. Conversely, assume there is 

a sequence {t }, t >0,t 0 with x = xt for each n 

n J n J n n 

and let t e R. If t = kt for some integers k and n, then 

n 

xt = x(kt ) = x by the lemma. Otherwise since t >0 for 
n ' J n 

all integers n, t/t e R and is located between two integers 

k and k +1. Thus k < t/t < k +1 and 
n n n n n 

k t < t < (k + l)t . Let 6 = min {(k + l)t - t, t - k t } 

nn n n n n * nn 

Since t -*■ 0, there exists an integer m > n such that 

t < 6/2. As above there exists integers k and k + 1 
m mm 

such that k t < t < (k + l)t . Since this interval has 
mm m m 

length (k + l)t - k t = t < 6/2, this interval must be 
m m m m m 

contained inside the previous interval since 

max {(k + l)t -t, t - k t } < 6/2. Thus we have 

m m mm 

kt <kt < t <"(k +l)t <(k +l)t. 

nnmm m m n n 

Continuing in this manner, we have constructed a sequence 

{k t } with the property that k t t . Now by the 
nn t' t' j nn 



continuity axiom, xCk^t^) •> xt and since x(k t ) = x for each 
n, by the lemma, we have x = xt . Thus x is critical. 

To prove the equivalence of statements 1. and 5. we 

first note that if x is critical, then 5. follows trivially. 

On the other hand, suppose x = x[a,b] for some a < b and 

let t e R. If t = kt^ for some integer k and t^ in [a,b], 

then xt = x(kt^) = x by the lemma. Otherwise, since [a,b] 

is closed, there is a sequence {t^} in [a,b] with t -* a. 

Consider the sequence {t n - a}. Clearly (t - a) > 0, 

(t - a) -*■ 0 and x = xa = x(-la) by the lemma. Since 

t n e [a,b], x = xt . Combining these last two statements 

we have x = xa = x(-a) = xt (-a) = x(t - a). Thus the 

n n 

hypotheses of statement 6. are satisfied and x is critical. // 

From this theorem we know that if we follow a point of 
any trajectory in a dynamical system and discover that that 
point is at rest for any period of time, then the point 
remains at rest indefinitely. 

Continuing to seek motivation from our examples, we 
again refer to Example 2.3- We see that if x is any point 
on the trajectory, x = X(X^,t) for some t e R but also 
x = X(Xp,t + 2nir) for n = ±1, ±2, ... Transfering this to 
our notation of a dynamical system, we have that 
x = Xgt = Xg(t + 2mr). As we would expect, any point with 
this property is called a periodic point. Formally, we 
state : 
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DEFINITION 2.17 A point x e X is said to be periodic if 
there is a T / 0 such that 

xt = x ( t + T) for all t e R. 



A number T for which this holds is called a period of x. 

If x is a periodic point then both the motion and the 
trajectory are said to be periodic. 

The restriction that T ^ 0 is necessary since it is 
true that xt = x(t + 0) for all t but the point x may not 
be periodic. It is also noted that any critical point x 
is periodic since every T e R is a period for x. The 
following theorems and remarks prove some of the expected 
properties of these points. 

THEOREM 2.1 8 If {x^} is_ a sequence of periodic points with 
positive periods T^ -*■ 0 and x n x, then x is_ critical . 

PROOF As in the proof of Theorem 2.16, given any t e R, 

there are integers k and k +1 such that 
& n n 



k T < t 
n n — 



< (k 



n 



+ 1 )T = 
n 



k T + T , 
n n n 



Since T ->■ 0, k T -*■ t . Then x = x T since x is periodic 
n’nn nnn n 

and x T = x (k T ) by lemma 2.15. Combining these last 
n n n n n J 

two statements we have x = x T = x (k T ) -*• x t. 

n n n nnn n 



Now 



// 



since x x we have x = xt. As t e R was arbitrary, 
x = xt for all t e R. Hence x is critical. 

From the definition we can see that a point x e X is 
periodic if and only if there is a T ^ 0 with x = xT. For 
if x is periodic, we know there exists a T f 0 such that 
xt = x(t + T) for all t e R. But in particular, for t = 0, 
we have x = xT. On the other hand, if there exists a T ^ 0 
such that x = xT and if t is any element in R, then 
xt = xT(t) = x(T + t) = x(t + T) . Since t was arbitrary, 
we see that x is periodic. 

In Example 2.3 we noted that T = 2-rr was a period of x 

and T = 4-rr was a period of x. Clearly in this example, 

x has an infinite number of periods. We would expect that 

x would have a smallest non-negative period, which we will 

call its fundamental period. If there exists an infinite 

number of periods for a point x in a dynamical system whose 

fundamental period is 0, then as an immediate consequence of 

the last theorem, .we know that x is critical. This is true 

since we can certainly form a sequence {T^} of the periods 

such that T ->-0 and since the sequence {x } where x = x 
n n n 

for all n has the property that x -»• x, the hypotheses of 
the theorem are fulfilled. The following theorem proves 
the existence of a fundamental period and also proves why 
we need only consider positive periods. 

THEOREM 2.19 If x e X is_ periodic but not critical , then 
there is a T > 0 such that T is_ the smallest positive period 
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of x. Further T' is_ any other period of x, then 
T* = nT for some integer n. 

PROOF Let P = {t > 0: t is a period of x}. P j- 0 since if 

T is a period of Xj T / 0, then x = xT. Thus 

x(-T) = xT(-T) - x(T-T) = x so -T is also a period of x. 

Since either T or -T is positive, P ? 0. Now set T = inf P. 

Such a T exists since P is bounded below by 0 and hence 

has an infimum. Clearly T > 0 so suppose T = 0. Then there 

is a sequence {t } in P with t 0. Since x = xt for 

all n, x must be a critical point by Theorem 2.19. Hence 

T > 0. We see that T is also a period of x since there is 

a sequence {t } in P with t -* T because T = inf P and 
n n 

x = xt for all n. Thus by the continuity axiom, 
x = xt -> xT. Thus T is a period of x and by definition 
it is the smallest positive period of x. Finally let t e R 
be any period of x. If t ^ nT for any integer n, then 
there exists an integer n such that nT < t < (n + 1)T. 

But since T is a period of x, nT is a period of x by 
lemma 2.15. Thus we have x = xt = x(nT) which gives 
xt(-nT) = x(nT)(-nT) = x(nT-nT) = x. Then 

x = xt(-nT) = x(t-nT) which shows that t-nT is also a period 

of x. But since 0 < t-nT < T we have a contradiction 

since T was the smallest positive period of x. // 

REMARK 2.20 If x e X is periodic, then y(x) and y + (x) 
are compact . 



PROOF Let x be periodic with positive period T. Then 
Y(x) = x[0 jT] and y + (x) = x[0,T] since for any t e [0,T], 
xt e x[0,T] and for any t with t > T, t = t^ + kT for some 
integer k and t^ e [0,T]. Clearly similar results hold for 
t < 0. From Remark 2.6, x[0,T] is compact. Thus y(x) and 
y + (x) are compact. // 

From this analysis it is now clear that if any tra,1 ectory 
in a dynamical system ever crosses itself , then the 
traj ectory must be periodic . For suppose xt^ = xt£ for 
t^ < t 2 « Then x = x(t 2 - t^) = xT where T = t 2 - t^ > 0. 

But this shows that x must be periodic. 
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III. PROLONGATIONS AND LIMIT SETS 



Referring again to Example 2.2, it was noted that in the 

damped spring case the solution converges to the point (q) 

as t -*■ «>. This solution converges in the sense that there 

is a sequence {t } with t e R for all n, t -*■ «> and 

n n n 

xt n (q). This point is called a positive limit point. 
Formally, the definition is stated: 

+ - X 

DEFINITION 3.1 Define maps A , A from X into 2 by setting 
for each x e X, 

1. the set A + (x) = {y e X: there is a sequence 

{t } in R with t ->-+«> and xt y} , 
n n n J 5 

2. the set A~(x) = {y e X: there is a sequence 

{t } in R with t -*■-«> and xt -»■ y}. 
n n 

For any x e X, the set A + (x) is called its positive limit 
set , and the set A - (x) is called its negative limit set . 

REMARK 3.2 If x e X is periodic, then A + (x) = A (x) = y(x). 

PROOF Let x e X be periodic. By Remark 2.20, y(x) is 

compact. Let y e y(x). Then y = xt for some t e R. Since 

x is periodic, y = x(t + T) = xT(t + T) = x(t + 2T). Clearly 

this idea can be extended inductively so consider the 

sequence {t + nT}. Since T > 0, t + nT ■+ +°° and y = x(t + nT) 

for all n. Thus x(t + nT) ■+■ y and y e A + (x). Hence 
+ + 

y(x) C. A (x). Now let y e A (x). Then there exists a 
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y. For all n, 



sequence {t } in R with t -+ -f 0 ° and xt 
n n n 

xt n e y(x). Hence {xt^} C. y(x). Since this sequence 
converges. It must converge to a point in y(x) because 
y(x) is compact. Thus y e y(x) and A + (x)C y(x). Hence 
A + (x) = y(x). The proof for A~(x) = y(x) is entirely 
analogous. // 

An immediate consequence of this remark is that i_f a 
point x is. a critical point , then A + (x) = A - (x) = x. From 
these two observations, it is an easy task to determine the 
positive and negative limit sets of Examples 2.1, 2.2, and 
2.3«, because in Example 2.1 we have critical points, 
periodic trajectories, and trajectories which clearly do 
not converge. In Example 2.2, we have either a periodic 
trajectory for the undamped case or a trajectory whose 
positive limit set is the origin and negative limit set is 
empty. Finally in Example 2 . 3 we have a periodic orbit. 
However, a point may not be periodic and yet have more than 
one point in its positive or negative limit set. The 
following examples help to illustrate this idea. 

EXAMPLE 3*2 Consider the differential system defined in 
2 

R by the following equations in polar coordinates: 



dr 

dt 



d6 

dt 



r (1 



1 . 



r) 
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An extensive analysis is not necessary here because all the 

facts needed can be derived directly from the equations. 

This system does in fact define a dynamical system [Bhatia, 

Szego, p. 20]. Since the change in 6 is constant, only the 

change in r affects the system. Clearly for r = 0, 
dr 



dt 



= 0 so the origin is a critical point. For r = 1, 



dr _ 



dt 



= 0 so r is again constant with time and the solution 



describes a periodic trajectory. Suppose r < 1. Then 
dr* 

is positive so r is increasing with time. If r > 1 , 



dr* 

then is negative and r is decreasing with time, 



Such 



trajectories are shown in the diagram below. It is noted 
dr 

that although ^ is always positive for r < 1, r will never 



reach r = 1, for as r -*• 1, 



dr 

dt 



0 . 




As this diagram shows, for points p with 0 < r < 1, A + (p) 
is the unit circle and A~(p) is the origin. For points p 
with r > 1, A + (p) is the unit circle and A (p) = 0. These 
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facts become clear when it is realized that a sequence 



{t^} can be found with t n -*■ «. such that xt^ approaches any 
given point on the unit circle. Since the unit circle is 
periodic and the origin is a rest point, their positive 
and negative limit sets are obvious. The direction of 
the trajectories is counter-clockwise because 0 is increasing 
with time. 

EXAMPLE 3.3 Consider the differential system defined in 
? 

R by the differential equations 



x = f (x,y ) 



y = g(x,y) 



where the functions f and g are defined by 



f(x,y) = 




0 if | x | > 1 



,y (1 - rh if | x | < 1 

(1 + y J (1 - p ( x)q (y ) ) 



g(x,y ) = -1 if x £ -1 

x if | x | < 1 . 




1 if x > 1 



Here the functions p(x) and q(x) are any continuously 
differentiable functions satisfying the following conditions: 



p(x) = 0 
P(x) = J 
0 < q (y ) < fjj- 
q(y) = 0 



if 0 < x < 1 , 
if x < 0, 
if x >_ 1 , 
if y < o, 
if y > 0. 



This system defines a dynamical system [Bhatia, Szego p. 21] 
and it is analyzed in the same manner as the preceding 
example. For x = y = 0, f(x,y) = 0 and g(x,y) = 0 so x and 
y are unchanging with time. Consequently, the point (q) is 
a critical point. For x _> 1, f(x,y) = 0 and g(x,y) = 1. 

The trajectories of such points are straight lines parallel 
to the y-axis with direction in the +y direction. A similar 
analysis shows that for x <_ -1, the trajectories are straight 
lines parallel to the y-axis with direction in the -y 
direction. Now suppose 0 < x < 1. Then g(x,y) > 0 so y 
is increasing with time. If y < 0, f(x,y) > 0 so x is 
increasing. If y > 0, then f(x,y) < 0 so x is decreasing. 

A similar analysis holds for -1 < x < 0. Thus the solutions 
describe spiralling trajectories. The question is do they 
spiral outward or toward the origin or are they periodic 
trajectories. Such an analysis from the equations is very 
involved but the trajectories do in fact spiral outward 
[Bhatia, Szego, p. 22]. The trajectories are sketched in 
the diagram below. 



5 1 * 



1 




Prom the diagram one can see that for any point except the 

origin in the strip -1 < x < 1 , its positive limit set is 

the two straight lines x = ±1 and its negative limit set 

is the origin. Since the origin is critical, its positive 

2 

limit set is itself and for any other point in R , its 
positive and negative limit sets are empty. 

Proceeding now with the concept of limit sets, some very 
important theorems are useful. 

THEOREM 3. 4 For^ an£ x e X, 

1. the sets A + (x) and A~(x ) are closed and invariant , 

2. the closure cl (y + (x)) = Y + (x) U A + (x) and the 
closure cl (y”(x)) = y~(x) U A - (x). 

PROOF 1. Consider the case of A + (x). Let {y^} be a 
sequence in A + (x) with y^ y. If it can be shown that 
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y e A + (x), then A + (x) is closed. Proceeding in this manner, 

it is noted that for each positive integer k, there is a 

sequence {t^} in R with t^ -*■ + 00 and xt^ -*• y. since each 
n n n J k 

y^ e A + (x). We may assume without loss of generality that 

k k 

d(y ,xt ) < 1/k and t > k for n > k where d is the metric 
k n n 

for the space X. These assumptions may be made since the 
xt get arbitrarily close to y. and hence any point in the 
sequence at distance greater than 1/k from y^ may be removed 
with the resultant subsequence converging to y^. The same 
reasoning applies to t £ k. Consider now the sequence 
{t^Mn R with t n = t^. The significance of this "diagonal" 
sequence is in consideration of (xt^} and for purposes of 
clarity, we sketch the latter sequence below:- 



V.. V V.. . 




Then t +°° and the claim is that xt n + y. To see this, 
observe that 

d(y,xt n ) £ d(y,y n ) + d(y n ,xt n ) < d(y,y n ) + 1/n. 

Since 1/n and d(y,y n ) tend to zero, this implies that 
d(y,xt ) -»■ 0. Consequently, xt n y and y e A + (x). Thus 
A + (x) is closed. To see that A + (x) is invariant, let 
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y e A (x) and t e R be arbitrary. There is a sequence 

{t^} in R with t^ +<» and xt -> y. Then by the continuity 

axiom, xt (t) yt. Since xt (t) = x(t + t) and 
n n n 

(t n + t) +oo this implies that yt e A + (x). Thus A + (x) is 
invariant . 

To prove the second part, it is recalled that y + (x) = xR + . 

By the definition of A + (x) it follows that 

cl (y + (x) ) 3 y + (x) V A + (x). Now let y e cl (y + (x)). Then 

there is a sequence {y n l in y + (x) such that y n y. Since 

this sequence is in y + (x) it follows that y^ = xt^ for some 

t e R. Either the sequence {t^} has the property that 

t •> +°°, in which case y eA + (x), or by the Bolzano-Weierstrass 

theorem there is a subsequence {t } with t t e R + . 

n k n k 

But then xt xt which is in y + (x), and since also 

n k 

xt -»• y, it follows that y = xt e y + (x). Thus 
n k 

+ + + 

cl (y (x))C y (x)U A (x). The proofs for the negative 
case are entirely similar. // 

The properties of limit sets just proved will be found 
to be very useful in further work. 

THEOREM 3-5 For x,y e X, if y e y + (x), then A + (y) = A + (x). 

PROOF Let y e y + (x). Then there is a t e R such that 

y = xt . Now let u e A + (y). Then there exists a sequence 

{t }in R with t •> +<» such that yt -* u. Consider the 
m m m 



57 













•I 






sequence {t + t } in R. Clearly t + t -*• +®. Then 

m m 

yt -*■ u Implies yt = xt(t ) = x(t + t ) -*• u. Thus 
m m m m 

u e A (x) and hence A + (y)C. A + (x). 

On the other hand, let z e A + (x). Then there is a 

sequence {t } In R with t +°° and xt z. Consider the 
m mm 

sequence (t - t) in R. Clearly (t - t) -*■ « since t -*■ «, 
ni m m 

Then xt -*• z implies xt = x(t + t - t) = xt(t - t) = 
n n m m 

y (t m - t) -*• z. Thus z e A + (y) and A + (x) C A + (y). 

REMARK 3.6 Let x,y e X such that x is not periodic. Then 
the set 0 {y + (y ) : y e y + (x)} = 0. 



+ + 

PROOF Since y e y (x), there is a t e R such that y = xt . 

+ + + 

Now suppose s e A(y (y): y e y (x)}. Then s c y (y) which 

Implies that s = yt^ for some t^ e R + . But since y = xt , 

it follows that s = xt(t^) = x(t + t^). Let tj e R such 

that tg > t + t^. Then for y = xt 2 , s l y + (y). Thus 
+ + 

s t fliy (y): y e y (x); and we have a contradiction. Hence 
fiW + (y) : y e y + (x)> = 0. // 

REMARK 3.7 For any x e X, A + (x) =fi(cl (y + (y)): y e y + (x)} 
my + (xn): n is an integer}. 



PROOF If x is periodic, the proof is trivial since 
cl (y + (y)) = A + (x) for all y. So we assume x is not periodic, 
Let {t^} be any sequence in R such that t n +». Without 
loss of generality it may be assumed that t > 0 for all n 
since any points in the sequence for which this does not 
hold may be removed from the sequence without affecting 
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t n ** + 00 * Since the intersection in the statement is taken 

over all possible y's in Y + (x), set y = xt where each 

n n 

t n e tt n ). Now Y + (y n ) = {y n t: t e R + } = {x(t n + t): t e R + }. 
The claim is that A {cl ( y + ( y > ) } = A (cl (Y + (y n )))- 

By Theorem 3-4, A (cl ( y + ( y ) ) > = A{y + (y)U A + (y)} = 
A{Y + (y))^ A { A + (y ) } and 

n{cl (Y + (y n ))> = n{ Y + (y n ) U A + (y n )} =n( Y + (y n )}^ A { A + (y n ) } . 
By Remark 3.6, A{Y + (yH = J3 = A’(Y + (y n )} and by Theorem 3*5, 
A + (y) = A + (x) = A + (y n ). Thus 
A {cl (Y + (y))) = A { A + (y ) } = A + (x) and 

A {cl (Y + (y n ) ) ) = A{A + (y n )} = A + (x). Then, as claimed, 

A {cl (y 4 (y))} = A {cl ( y + ( y n ) ) ) and as shown above, 

A {cl ( y + ( y ) ) : y c y + (x)} = A + (x). To prove the second 
part, let the sequence {t^} = {n}. Then the proof is 
identical. // 

In the following proof and elsewhere in this thesis, 
the notation S(M,e) and S[M,e] will be used for the sets 
{x: d(x,M) < e}, and {x: d(x,M) <_ e} respectively where 
M C x, e e R and d is the metric for the space X. In the 
case where M = {x}, the notation will be shortened to 
S(x,e) and S[x,e]. 

THEOREM 3.8 Let x e X. If cl (y + (x)) is compact , then 
A + (x) is_ a non - empty compact and connected set . 

PROOF Let cl (y + (x)) be compact and let {t- n > be any 
sequence in R + such that t -*• +°°. Now consider the sequence 
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{xt n >. Since xt n e cl (y + (x)) for all n and cl (y + (x)) 
is compact, this sequence has a convergent subsequence 
{xt k >which converges to a point y e cl (y + (x)). Thus since 
t k +oo and xt k y, we have A + (x) ¥■ 0 since y e A + (x). 

Now cl (y + (x)) = y + (x) U A + (x). By Theorem A + (x) is 

a closed set and since it is a closed subset of a compact 
set, A + (x) is itself compact. 

The proof of connectedness is by contradiction. Suppose 
A + (x) is compact but not connected. Then A + (x) = P U Q 
where P,Q are non-empty, closed disjoint sets. Since 
cl (y + (x)) is compact, it is locally compact and hence there 
exists an e > 0 such that S[P,e] and S[Q,e] are compact and 
disjoint. To see that S[P,e] is compact, it is noted that 
in a locally compact space, every point has at least one 
compact neighborhood. Taking the union of all these compact 
neighborhoods, yields a compact neighborhood of the set P. 
Then there is an e > 0 such that S[P,e] is a closed subset 
of this neighborhood and hence S[P,e] is compact. The 
same reasoning applies to S[Q,e]. Since a metric space 
is normal, there exist disjoint open sets containing P and 
Q respectively. Hence for a proper choice of e, the sets 
S[P,e] and S[Q,e] are disjoint and compact. 

Now let y e P and z e Q. Then there exist sequences 

{t } and {t } with t -»-+<» and t +°° such that xt y 
n n n n n 

and xx n z. Without loss of generality it may be assumed 
that xt n e S(P,e), xt^ e S(Q,e) and x n - t n > 0 for all n. 
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Since the trajectory segments x[t ,t ] , n = 1 , 2 , . . . are 

n n 5 

compact connected sets, they clearly intersect S(P,e) and 
S(Q,e). Thus in particular, there is a sequence {T^} 
with t n < T n < t n such that xT n e H(P,e) where 
H(P,e) = {x: d(x,P) = e}. Since H(P,e) is a closed subset 
of S[P,e], H(P,e) is compact. Since {xT^} is a sequence 
in H(P,e), it has a convergent subsequence {xT m > which 
converges to a point y' and because +», it follows 

that y 1 e A + (x). But y' £ P U Q which contradicts the 
assumption that A + (x) = P U Q. Hence A + (x) is connected. // 

Since only the local compactness of cl (y + (x)) was used 
in the proof, an immediate corollary to this theorem is 
that if^ the space X is_ locally compact , then A + (x) is_ connected 
whenever it is compact . 

The converse of the previous theorem holds if the space 
X is locally compact. However, the proof relies on the one- 
point compact ification of X. The fact that the one-point 
compactif ication is an extension of a dynamical system to 
another dynamical system is proved in the following lemma. 

LEMMA Given a dynamical system (X,R,tv) on a locally compact 
space X, let X 1 = X U {w} be the one-point compactif ication 
of X and define tt ' : X'xR X* by 
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Then (X’jRjTt') is a dynamical system on X' . 



PROOF First, by definition of tt 1 



TT ' ( x , 0 ) = 




tt(x,0) = x if x e X 



if x = w . 



Thus tt' ( x,0) = x for all x e X*. 

Secondly, if x e X, then tt ' (it ' (x,t 1 )',t 2 ) = Tr’(x,t, + t 2 ) 



follows trivially so suppose x = to . Then 

tt' (tt* (co,t 1 ) ,t 2 ) = ir'(w,t 2 ) = to = TT'(cj,t 1 + t 2 ). Thus 

tt* ( tt 1 (x,^) ,t 2 ) = TT'(x,t 1 + t 2 ) for all x e X' . It remains 

to show that tt* is continuous. To see this, we will utilize 

the characterization that if the inverse image of an open 

set is open, then the map is continuous. 

So let U be any open set in X'. If to £ U, then U is 
open in X and thus (it ' (U) = tt”^(u). Now it”^(U) is open 

because tt is continuous and hence (tt')~^(U) is open. If 
to e U, then X-U is compact and since to £ X'-U, X'-U = X-U 
and hence (tt ' ) -1 (X' -U) = (tt') _1 (X-U) = tt _1 (X-U). Now since 
tt is continuous, tt“^(X-U) is closed and hence (tt* ) - ^(X'-U) 
is closed so (tt') _ ^(U) is open. Thus the inverse image of 
every open set is open and hence tt* is continuous. // 

THEOREM 3*9 If the space X is_ locally compact , then for 
any x e X, cl(y + (x)) is_ compact whenever A + (x) is_ non - empty 
and compact . 
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PROOF Let A + (x) be non-empty and compact. Consider the 
extension of the dynamical system to the one-point compacti- 
fication X' = X U {to}. Certainly, for any x e X, 

A (x) = A (x) U {to}or A (x) = A + (x) where A + (x) is the 
positive limit set of x in X'. However since A + (x) is 
compact, it follows that A + (x) = A + (x) holds because 

-f. -f- -L 

A (x)=cl^, (A (x)) = cl^(A (x)) = A (x). Also we have that 
for any x e X, y + (x) = y + (x) since 7r'(x,t) = ir(x,t). Thus 
we have 

cl x ,(-y + (x)) = y + (x) U A + (x) = y + (x) U A + (x) = cl x (y + (x)) 

and since cl x ,(y + (x)) is a closed subset of a compact space, 
it is compact. Thus cl(y + (x)) is compact. // 

THEOREM 3*10 If cl(y + (x)) is_ compact , then d(xt,A + (x)) -*■ 0 
as t -*• +<». 

PROOF Let cl(y + (x)) be compact and suppose d(xt,A + (x)) / 0 

as t + +oo. Then there exists a sequence {t^} in R such that 

t -»■ +oo and d(xt ,A + (x)) > e for all n > N where N is an 
n n J 

integer and e > 0. But since cl (y (x)) is compact, the 
sequence (xt n > has a convergent subsequence {xt^} which 
converges to a point p e cl (y + (x)). Now since t -*■ + 00 , 
it follows that t^ + 00 and since xt^-»- p, we have that 
p e A + (x). Also since xt^ + p, this sequence gets arbitrarily 
close to p and hence d(xt ,p) < e for k sufficiently large. 
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Hence we have d(xt k , A + (x)) < e for k sufficiently large 
but this is a contradiction and the theorem is proved. // 

The study of limit sets can be very useful in predicting 
the behavior of a system. If the system describes a physical 
problem and the trajectories are sketched, it is a relatively 
simple task to discover the limit sets and hence predict 
the future behavior of the system. For example, consider 
again the pendulum problem, Example 2.1. If it is known 
that the pendulum is started at the equilibrium position and 
given an initial angular velocity of to = 2k, then one would 
predict that the pendulum would approach the position where 
0 = it and would never change direction to swing back. This 
result would be known without having someone watch the 
pendulum to actually see what happens. It is true that 
this is a very elementary problem, but the ideas involved 
have broader applications to problems of much greater 
difficulty . 

It is important that the reader have a firm understanding 
and familiarity with the symbols and sets thus far defined 
because the following additional sets and symbols become 
quite important in the next two sections of this thesis. 

DEFINITION 3.11 For each x e X, define the following 
subsets D + , D _ , J + , J~ of X, 

1. the set D + (x) = {y e X: there is a sequence {x^} 

in X and a sequence {t } in R such that x x and x t -*■ y } , 

n n n n 
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2. the set D (x) = (y e X: there is a sequence (x^J 

in X and a sequence (t } in R - such that x + x and x t y}, 

n n n n 

3. the set J (x) = {y e X: there is a sequence {x^J 

in X and a sequence {t Jin R + such that x -*■ x, t +°°, 

and x n t n + yj, 

4. the set J (x) = {y e X: there is a sequence {x^J 

in X and a sequence {t Jin R~ such that x -*■ x, t 

n n * n * 

and x t yj. 
n n 

For any x e X, the set D + (x) is the first positive prolongation 
of x and the set D - (x) is the first negative prolongation 
of x. 

The sets J + (x) and J~(x) are called, respectively, the 
first positive and the first negative prolongational limit 
set of x. 

These sets are called "first" because higher-order 
prolongations can be defined. However, in this thesis, 
such higher-order prolongations will not be discussed. 
Consequently, the word "first" will be suppressed. 

From the definitions it is clear that any x e X, 

D + (x) 3 y + (x), D~(x) y (x), J + (x) O A + (x), and J (x)^ A (x). 
To achieve an intuitive feeling for these sets and to see 
that the inclusions just noted may be proper, we present 
the following example. 

EXAMPLE 3*12 Consider a predator-prey problem where it 
is assumed that a single predator population feeds on a 
single prey population. Let x denote the size of the prey 
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population and let y denote the size of the predator 
population. Then this problem can be expressed by the 
system 



dx 

dt 



= ax - bxy, x >_ 0, 



|| = cxy - dy, y > 0, 



where a,b,c, and d are all non-negative constants. 

These equations support our intuition since we would 

expect that any change in the size of the predator population 

depends not only on the size of the prey population but also 

on the size of the predator population. For example if the 

prey population is small and the predator population is 

large, we would expect that the predator population would 

decrease. This expectation is confirmed by the equations. 

Rather than solving this system in terms of t, it will 

be analyzed directly from the differential equations. 

Clearly for x=0=y, ^-=0=^- and the point (q) is 

a critical point. One can also see from the equations 

d / c 

that the point ( a / b ) is a critical point. For any point 
P = (n)> we have ^ > 0 and = 0. Thus its associated 
trajectory lies on the x-axis. Physically, this is the 
case when there is no predator population so the prey 
population increases. For any point p = (^), = 0 and 

< 0 so its trajectory lies on the y-axis with direction 
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in the negative y-direction. This is the case when there is 

no prey population so the size of the predator population 

decreases. For small x and large y, ^ < 0 and ^ < 0. 

Thus both x and y are decreasing and x continues to do so 

d x 

until y = a/b. At this point = 0 and as y continues to 
decrease, x begins to increase. When x reaches a value of 
d/c, = 0 and as x continues to increase, y begins to 
increase. The physical situation described by this case 
is quite clear. The trajectories for this system are 
sketched below. 



r 




For any point p = (^) we have A + (p) = (5) but 



j + (p) = {(*):x > 0}. The fact that A + (p) = (jj) means that 
if there is no prey population, the predator population will 
eventually decrease to nothing. The fact that J (p) = (q) 
means that even though the prey population may decrease 
considerably so that the predator population decreases to 
nothing, there will be enough of the prey population left 
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to eventually Increase Its size to any amount. We also 
notice from this diagram that for any point p = (°), 

A + (p) = (q) but J + (p) = { (q) : x > 0}. Also 

D + (p) = Y + (p) b {(q)}, D (p) = y~ (p ) and J~(p) = 0. For 

any point p = (£), D + (p) = y + (p), J + (p) = 0, D~(p) = y”(p)U 

{(y)} 5 and J (p) = {(^)). For any other point p 3 

D + (p) = Y + (p), D (p) = y (p), J + (p) = 0 = J _ (p). It may 

be argued here that we cannot form infinite sequences since 

"size" is not a continuous function. However, in an actual 

experiment, the following values for the constants a,b,c,d 

Q 

were observed: a ~ 1 , b - .01, c ~ 3 x 10 , and d ~ 0.45. 

These values were obtained by observing paramecia feeding 

d x 

on yeast cells. The value of y for which changes sign 
is then at y - 100 and the value of x for which changes 
sign is at x ~ 1.5 x 10^. These values are large enough 
to assume continuity. Also any mathematical model idealizes 
the physical situation: the question is how "close" [Bailey, 

p. 245, 246]. 

Now we consider some basic properties and important 
theorems about these limit sets. 

THEOREM 3-13 For an^_ x e X, 

1. the set D + (x) is^ closed and positively invariant , 

2. the set J + (x) is_ closed and invariant , 

3. the equality D + (x) = y + (x)U J + (x) holds . 
Analogous results hold for D~(x) and J~(x). 
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PROOF Statement 1 will follow from statements 2 and 3 so 
statement 2 is proved first. This proof parallels the proof 
of Theorem 3.^. Let {y^} be a sequence in J + (x) with 
y n y* For each integer k, there are sequences {x k } in X 
and{tjj ) in R + with x k + x, t k -*■ +°°, and x k t k -*■ y R . Without 
loss of generality, it may be assumed that t k > k, 

d (x k ,x) 1 1/k > and d(x n t n>yk ) - 1/k for n 1 k - Now consider 
the sequences {xJJ}, {tJJ}. Clearly x” + x, t” -»■ +°°, and 
x£t£ y. This is true since 

d(x n t n ,y ^ - d(x n t n ,y n ) + d(y n ,y) - 1/n + Thus 

+ + 

y e J (x) and J (x) is closed. 

To see that J + (x) is invariant, let y e J + (x) and t e R. 

There is a sequence (x n ) in X and a sequence (t^) in R + 

such that t -*■ + 00 , x -*• x, and x t -»■ y. Now consider the 

sequence {t + t). Clearly t + t -*• + 00 , and ' 

x(t + t) = xt(t)-»- yt. Since x + x, it follows that 
n n n n ' J n * 

+ 

yt e J (x). Since, t e R was arbitrary, it follows that 
J + (x) is invariant. 

Consider statement 3. Observe that 
D + (x) y + (x)Vj J + ( x ) always holds from the definition. So 
let y e D + (x). Then there is a sequence {x^} in X and a 
sequence {t } in R + with x -»• x, x t -*■ y. It may be assumed 
that either t t e R + or t n -*■ + 00 , if necessary by taking 
subsequences. In the first case x^t^ xt by the continuity 

axiom. Hence, xt = y e y + (x) as t e R + . In the second case 
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y e J + (x) by definition. Thus y e y + (x) ^ J + (x) and 
D + (x)C. y + (x) U J + (x). Therefore, D + (x) = y + (x) \J J + (x). 

To see statement 1, recall that cl (y + (x)) = y^(x) \J A + (x). 
But it has been shown that A + (x) C J + (x). Hence 
D + (x) = cl (y + (x)) U J + (x), because J + (x) "D A + (x). Since 
both these sets are closed and positively invariant, it 
follows that D + (x) is closed and positively invariant. // 

The following theorem will also be found useful in 
subsequent sections. 

THEOREM 3.1*1 Let X be locally compact . Then A + (x) ¥■ 0 

4 

whenever J (x) is non - empty and compact . 

PROOF Suppose A + (x) = 0. Then 

cl (y + (x)) = y + (x) U A + (x) = y + (x) since A + (x) = 0. Hence 
y + (x) is closed. It is claimed that y + (x) is disjoint from 

J + (x) for if y + (x) A J + (x) 0 0, then by the invariance of 

4* + *4* -f* 

J (x), y (x) C J (x). Since J (x) is compact, for any 

sequence { t } in R + such that t +<» } the sequence {xt^} 

must have a convergent subsequence. It follows that the 

point to which that subsequence converges is a point in 

A + (x) and hence A + (x) ^ 0. Thus y + (x) A J + (x) = 0 as 

claimed. Since J + (x) is non-empty and compact, 

d(y + (x) , J + (x) ) > 0. Thus since X is locally compact, there 

is a 6 > 0 such that S[J + (x),6] is compact and disjoint 

from y (x). Now choose any y e J (x). There is a sequence 

{x } in X and a sequence (t } in R + such that x x, 
n n n 
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t n +0o j and x^t^ -*• y . It may be assumed that 
x 0 S[J (x),S], x n ^ n eS[J(x),6] for all n. Then the tra- 
jectory segments x^O,^] intersect H(J + (x),6) v/here 
H(J (x),6) = {z e X:d(J + (x),z) = 6}. Therefore there is 
a sequence {T n >, 0 < x n < t n , such that x^r n e H(J + (x),6). 
Since H(J + (x),6) is compact, it may be assumed that 
x n x n -*• z e H(J + (x),6). By taking subsequences if necessary, 

it may be assumed that either t -*■ t e R + or t -*■ +°°. If 

n n 

T n -*• t e R , then by the continuity axiom ->■ xt = z 

so z e y + ( x ) which contradicts y + (x) f\ S[J + (x),6] = 0. 

If -> + 00 , then z e J + (x) , but this contradicts 

z e H(J + (x),<5) as J + (x) D H(J + (x),6) = 0. Thus the original 

assumption that A + (x) = 0 is untenable and hence 

A + (x) j- 0. // 

THEOREM 3*15 Let X be locally compact . Then J + (x) is 
non - empty and compact if and only if D + (x) is_ compact . 

PROOF Let J + (x) be non-empty and compact. Then by the 
preceding theorem, A + (x) is non-empty and, since it is a 
closed subset of J + (x), it is compact. But then cl (y + (x)) 
is compact by Theorem 3*9. Hence 

D + (x) = Y + (x) U J + (x) = cl (y + (x))\J J + (x) so that D + (x) 
is compact. On the other hand, suppose D + (x) is compact 

and assume J (x) = 0. But J (x) = 0 implies y (x) is 

+ + + 

compact. Then y (x) = cl (y (x)) so cl (y (x)) is compact 
but this implies A + (x) 0 0 by Theorem 3.8 and the assumption 
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that J + (x) = 0 is contradicted. Thus J + (x) 0 0 and since 
J + (x) is a closed subset of D + (x), it follows that J + (x) 
is compact. // 

Before proceeding to the next theorem, the set D + (M) 
must be defined where M <£. X. Thus D + (M) = {y e X: there is 
a sequence {x^} in X and a sequence { t n } in R + such that 
x x e M, and x t + y } . 

THEOREM 3 • 1 6 For any compact set M C X , the set D + ( M ) is 
closed and positively invariant . 

PROOF Let {y^} be any sequence in D + (M) such that y n + y. 

k k 

Then for each y^ there exist sequences {x^} in X and {t } 

+ k k k 

in R such that x v x, and x t ->■ y, where x. e M. Without 

n k n n k 

k k 

loss of generality it may be assumed that d(x n t n ,y k ) < 1/k. 

Now consider the sequences {x^} in M and { t in R + . Then 

x n t n y and since M is compact, it may be assumed that 
n n J n 

x n x e M. Thus x n t n + y since 
n n n J 

d(xJJtJJ,y) <_ d(xJJtJJ,y n ) + d(y n ,y) < 1/n + d(y n ,y). Then 

by definition, y e D + (M) and hence D + (M) is closed. To 

see that D (M) is positively invariant let y e D (M) and 

let t e R be arbitrary. Now y e D (M) implies there exist 

sequences {x } in X and {t } in R + such that x x e M 
n n n 

and x t + y. Now consider the sequence {t + t} in R . 
n n n 

Since x t -*■ y, it follows that x t (t) = x (t + t ) -^ yt. 
nn nn nn 

Thus yt e D + (M). Since t e R was arbitrary it follows that 
D + (M) is invariant. // 
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We conclude this section with an important theorem. 



THEOREM 3.17 Let x,y in X. Then y e J + (x) if and only 
if x e J”(y ) . 

PROOF Let y e J + (x). Then there is a sequence {x^} in X 

and a sequence {t } in R such that x -*• x, t +°°, and 

^ n n ’ n ’ 

x t y. Set x t = y , and x = -t . Then y^ -*• y, 
n n n n n’ n n n 



x -co and yx = xt(x) = x(t - t ) = x^ 
n J nn nn n nn n n 



x. 



Consequently x e J~(y). The converse holds in the same 
manner. // 
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IV. RECURSIVE AND DISPERSIVE CONCEPTS 



Using the sets defined in the last section, we continue 
with our analysis of dynamical systems by studying recursive 
and dispersive concepts. To motivate the idea of recursive- 
ness, consider again Example 3-3. Suppose we consider a 
set A to be the union of two circles centered respectively 
at the points (~q) and (q) with radius r where 0 < r < 1. 

In this example, all trajectories in the strip -1 < x < 1 
spiral outward from the origin with the lines x = -1 and 
x = 1 as their positive limit set. Thus for any point in 
the strip -1 < x < 1, its trajectory will eventually 
intersect A. The trajectory will pass out of A again but 
it nevertheless will return some time later and repeat the 
process. This is shown in the following diagram. 




The set A is said to be positively recursive with 
respect to the trajectory. We now state the formal definition 
of recursiveness. 
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DEFINITION 4.1 A set A C X is said to be positively 
recursive with respect to a set BC X if for each T e R 



there is a t > T and an x e B such that xt e A. Negative 
recursiveness may be defined by using the inequality t < T. 

A set A is self positively recursive whenever it is positively 
recursive with respect to itself. 

The simplest example of a self recursive set is the 
set {x} where x is a periodic point. By placing restrictions 
on the system by demanding that the entering of points from 
B into the set A happen with some regularity, other concepts 
such as almost periodicity can be discussed although such 
is not our purpose here. We now define some concepts and 
establish properties in connection with the concept of 
recursiveness . 

DEFINITION 4.2 A point x e X is said to be positively 
Poisson stable if every neighborhood of x is positively 
recursive with respect to {x}. 

This definition means that the trajectory through 
the point x eventually intersects every neighborhood of x. 

The following diagram illustrates such a situation. 
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To give a better understanding and some characteriza- 
tions of positively Poisson stable points, the following 
theorem is presented. 

THEOREM 4.3 Let x e X. Then the following are equivalent . 

1 . The point x is_ positively Poisson stable . 

2 . Given a neighborhood U of_ x and a T > 0 , 
we have xt e U for some t > T, 

3. The point x i_s contained in it s positive limit 
set A + (x) , 

4. The set cl (y + (x)) = A + (x), 

5. The set y(x) C. A + (x), 

6. For every e > 0 there is a t >_ 1 such that 
xt e S(x,e). 

PROOF The equivalence of statements 1- and 2 follows from 

the definition and our previous discussion. To see the 

equivalence of statements 2 and 3> consider a s'equence 

{T } in R + such that T +°°. To each n associate the 
n n 

neighborhood S(x ,1/n) of x. Then for each t there 

exists at > T such that xt e S(x,l/n). Clearly 
n n n 

t -*■ +~ and xt x. Thus x e A + (x). Conversely, let U be 
n n 

a neighborhood of x. Since x e A + (x) there is a sequence 

{t } in R + with t +~ such that xt -»■ x. Thus xt must 
n n n n 

eventually be in every neighborhood of x so that, in 

particular, given any T > 0 there exists a t^ > T such that 

xt e U. Next * we establish the equivalence of statements 
n 
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3 and 4. Let x e A + (x). Then since A + (x) is positively 
invariant, we have y + (x)C A + (x). Since 

cl (y (x)) = Y + (x) U A + (x), it follows that cl (y + (x)) = A + (x). 

Conversely, let cl (y + (x)) = A + (x). But since we know that 

cl (y + (x)) = Y + (x) U A + (x), it follows that Y + (x)d A + (x) 

and hence that x e A + (x). To see the equivalence of 

statements 3 and 5, let x e A + (x). Then, since A + (x) is 

invariant, xt e A + (x) for all t e R. Thus y(x) C A + (x). 

Conversely, let y(x) cz A + (x). Then for any t e R, 

xt e A + (x). In particular, for t = 0 we have x e A + (x). 

Finally, we prove the equivalence of statements 3 and 6. 

To this end let x e A + (x) and let e > 0. Then there exists 

a sequence (t } in R + with t -*■ 00 such that xt x. 

n n n 

Without loss of generality, we may assume t >_ 1 for all n. 

Then since xt x, this sequence must eventually be in 

every neighborhood of x. Hence there exists a t fi > 1 

such that xt n e S(x,e ). Conversely, choose a positive 

null sequence {e } such that e 0. Then for each n, we 

have a t >_ 1 such that xt n e S(x,e n ) . Since e n + 0, 

we have that xt x. Consider the sequence (t^} in R + . 

n n 

Either t -»■ +" or it contains a convergent subsequence 

(t } such that t -*• t e R + . If t ^ + 00 , x e A + (x) by 
n k n k n 

definition. If t t , then xt ■+ xt . But since xt x, 

n k n k n 

we have x = xt and therefore x is periodic with period t. 

It follows that x s A + (x) and the proof is complete. // 
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Prom this theorem it should be clear that whenever a 
point x is positively Poisson stable then so is xt for every 
t e R. This is true because if x is positively Poisson 
stable, then x e A + (x) which implies xt e A + (x) for all 
t e R by the invariance ofA + (x). But A + (x) = A + (xt) for 
all t e R by Theorem 3*5. Hence xt e A + (xt) and xt is 
positively Poisson stable. 

According to the above theorem, a point x is positively 
Poisson stable if x e A + (x), Similarly, a point x is 
defined to be negatively Poisson stable if x e A - (x). It 
is said to be Poisson stable if it is both positively and 
negatively Poisson stable. If a point x is Poisson stable, 
then both its motion and its trajectory are said to be 
Poisson stable. 

Prom the preceding theorem observe that if 
cl (y (x)) = A (x), then x is positively Poisson stable. 

Now consider the situation where y + (x) = A + (x). In that 
case the point x is positively Poisson stable since 
cl (y + (x)) = y + (x) \J A + (x) = A + (x). We will now establish 
that y + (x) = A + (x) if, and only if x ±s_ a periodic point . 

For if x is a periodic point, it is clear that y + (x) = A + (x). 
On the other hand, if y (x) = A (x), then x e A (x) and 
by the invariance of A + (x), y(x) = A + (x). Then 
xt * e y + (x) for t' < 0 and therefore there exists a T >_ 0 
such that xt ' = xT. Hence xt’(t) = xt(t') = xT(t). 
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Therefore xt = x(t + T - t') for all t e R, showing that 
x is periodic with period T - t' >0. 

With this last fact it is a natural question to ask 
whether there are Poisson stable points which are not 
periodic. The following example shows that such points do 
exist . 

EXAMPLE 4.3 Consider the dynamical system defined on a 
torus by the system 

dx w v 

at 

= af(x,y) , 

where f(x,y) is periodic in both x and y with period 1. 

That is f (x,y ) = f(x + l,y + 1) = f(x + l,y) = f(x,y + 1). 

Further, we have that f(x,y) > 0 for x and y not both 0 mod 1 

and f(0,0) = 0. The constant a is an irrational number. 

Clearly the point p = (q) is a critical point because 

at this point = 0 = Solving this system simultaneously 

d x 

we have — = 1/a which implies y = ax + c where c is some 
constant of integration. Thus the trajectories describe 
straight lines of slope a in the plane. We may assume 
a > 0. Nov; no trajectory can contain p except y(p) = (p) 
since p is a rest point and since ^ and ^ are both positive 
for all other points. 



79 



We note now that the function y = ax + c is not defined 
at the point p since at that point So we claim 

that there is only one trajectory which approaches the 
point p as t + The following diagram will help us to 

determine this trajectory. For purposes of this diagram, 
we assume 0 < a < 1. 




Consider the trajectory y = ax + (1 - a). As x + 1, it 
is clear that y 1. Hence this trajectory approaches 
the point (^) as t ->• °°. But on the torus, the point 
p = C° 0 ) = Hence this trajectory approaches p as t •* °°. 

Generalizing this approach to the case of an arbitrary 
positive value for a, the trajectory which approaches p as 
t oo is given by the equation y = ax + (1 - (a mod 1)). 

To verify this fact, let x -*• 1. Then y -*■ a + (1 - (a mod 1)) 
= 1 + (a - (a mod 1)). But a - (a mod 1) is just some 
integer n and hence as x + 1, we have y + n + 1. Thus 
as t this trajectory approaches the point ( n 1 ) but 

again, in the torus, this point and the point p are identi- 
fied. Thus this is precisely the trajectory which approaches 
p as t + “. Call this trajectory y 2 * Similarly there is 
only one trajectory which approaches the point p as 
t It is given precisely by y = ax + 0. It is clear 
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that as x 0, we have y •+• 0 and this verifies that this 
trajectory does in fact approach p as t -*■ -°°. Call this 
trajectory y^. 

Now we claim that no trajectory in the torus is periodic. 
If there is such a trajectory, then it will pass through 
the same point twice and without loss of generality, we 
may assume that this point lies on the y-axis. VJe may 
make this assumption because the trajectory must intersect 
the y-axis. Assume that this point is the point (®) where 
0 < c < 1. Then if the trajectory given by y = ax + c 
intersects this point twice, we have ( c ) = ( &n + c )mod 1 

for some n > 0. But this implies that c = (an + c) mod 1. 

Since 0 < c < 1, we know that c = c mod 1 = (m + c)mod 1 

where m is any integer. Hence it follows that 

(an + e)mod 1 = (m + c)mod 1. This implies that an = k for 
some integer k. Now since n / 0 , it follows that a = k/n. 
But this implies that a is a rational number which is a 
contradiction. Thus no trajectory can pass through the 
same point twice and hence no point except p is periodic. 

We also note here that the trajectories Y 1 and are 
not the same for if they were, then the equations 
y = ax + (1 — ( a mod 1)) and y = ax + 0 must be equal. 

But this implies that a mod 1=1 and thus a must be an 
integer which contradicts the fact that a is irrational. 

We now sketch several of the trajectories in the 
following diagram. 
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To complete our analysis of this example, we see that 
A + (y£) = {p} and that A (y-^) = (p). For the point p, we 
have A + (p) = A (p) = { p } . For any other trajectory, we 
have A + (y) = A (y) = the torus since these trajectories 
are not periodic. Also, A + (y 1 ) = A~(y 2 ) = the torus. 

Hence we have that all the points on y^ are positively 
Poisson stable since for each x e y^, x e A + (x). Also 
points on are negatively Poisson stable. All other 
points are Poisson stable. 

If we change this example so that f(x,y) > 0 for all 
x and y, then we will have no rest points and no periodic 
trajectories. Hence every trajectory is dense in the torus 
and the positive and negative limit sets of each point is 
the torus. 

We now introduce the notion of a non-wandering point. 

DEFINITION 4.4 A point x e X is said to be non - wandering 
if every neighborhood U of x is self positively recursive. 
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It is clear from this definition that if a point is 
positively Poisson stable then it is non-wandering. This 
idea can also 'be seen in the following theorem which 
gives some characterizations of non-wandering points. 

THEOREM 4.5 For any x e X, the following are equivalent . 

1. The point x is_ non-wandering , 

2 . The point x is_ ari_ element of its positive 
prolongational limit set J + (x), 

3. Every neighborhood of x is_ self negatively 
recursive , 

4. The point x is_ an element of its negative 
prolongational limit set J~(x). 

PROOF Assume x is non-wandering. Consider a null sequence 

{e }, 0 < e , e 0, and a sequence{t } in R with t -*• +°°. 
n n* n ’ n n 

Since each S(x,e n ) is self positively recursive, we have an 

x e S ( x, e ) and ax > t with x x e S(x,e ). Since 
n n nn nn n 

e + 0 we have x x and x t x and since t -► +°° we 
n n * n n n 

conclude that x e J + (x). Thus statement 2 holds. Assume 

X e J + (x). Then there is a sequence in X and a 

sequence {t } in R with x -*• x and t -► + 00 such that 
M n n n 

x t -*■ x. Now for any neighborhood U of x and T > 0 
n n 

there is an N such that t n > T, x n e U, and x R t n e U for 
n >_ N. Thus U is self positively recursive and, consequently, 
x is non-wandering. The equivalence of statements 3 and 4 
is proved in the same manner. It remains to show that 
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statements 2 and 4 are equivalent. This follows directly 
from Theorem 3.17* // 



We shall now prove some theorems which establish the 
connection between Poisson stable points and non-wandering 
points . 

THEOREM 4 . 6 Let x e X. Every y e A + (x) is_ non - wandering . 

PROOF It is necessary to show that if y e A + (x) for some 

x e X, then y e J + (y). If y e A + (x), then there is a 

sequence {t } with t -»- +» and xt -*■ y. Since t ■+■ +°° 
n n n n 

we may assume, if necessary by taking a subsequence, that 
t n+ ^ - t >_ n for each n. Then setting 

t = t ,, - t and xt = x we have x y, x(t,,-t) = 

n n+1 n n n n J * n n+1 n 

xt n+1 -*• y, and T n = t R+1 - t R -> +~. Thus y e J + (y) and 
y is non-wandering. // 

THEOREM 4 . 7 Let PCX be such that every x e P i_s either 
positively or negatively Poisson stable . Then every 
x e cl P is_ non - wandering . 

PROOF Let {x } in P and x + x. We must prove that 
n n 

x e J + (x). For each n Me have that either x n e A + ( x n ) or 

x^ e A”(x ). Thus by taking a subsequence if necessary, we 

may a'ssume that x e A + (x ) for all n or that x^ e A (x ) 

n n n n 

for all n. In the first case, for each n there is a 

t > n with d(x ,x t ) < 1/n. Then clearly 
n n n n 

d(x, x nt n ) 1 d(x,x n ) + d ( x n » x n t n ) £ d ( x > x n ) + l/ n • Th i s 
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shows that x n t^ -*• x and consequently x e J + (x). In the 
second case similar considerations show that x e J“(x). 

Thus by Theorem 4.5, every x e cl P is non-wandering. // 

A partial converse of this last theorem will now be 
presented without proof. The proof is entirely constructive 
but not particularly instructive. A proof may be found 
in Bhatia and Szego. 

THEOREM 4 . 8 Let X be_ complete and let every x e X be 
non-wandering . Then the set of Poisson stable points P 
is dense in X. 

To see how several of these concepts are related, we 
present the following summary. It is clear that a critical 
point is a periodic point. Also, all periodic points are 
Poisson stable since the positive and negative limit sets 
of a periodic point are the periodic trajectory. Certainly 
the periodic point is contained in its periodic trajectory. 
Finally, all Poisson stable points are non-wandering points 
points for if x is a Poisson stable point, then 
x e A (x)C J (x). These ideas are summarized in the 
following diagram. 
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We now define a different type of stability. This 
notion will become important in discussing "dispersive" 
concepts . 

DEFINITION lJ.9 For any x e X, the motion it is said to 

x 

be positively Lagrange stable if cl (y + (x)) is compact. 
Further, if cl (y (x)) is compact, the motion tt is called 

X 

negatively Lagrange stable . It is said to be Lagrange stable 
if cl (y(x)) is compact. 

If X = R n , then the above statements are equivalent to 
the condition of the sets y + (x), y”(x), y(x) being bounded, 
respectively. The following statements about this notion 
are easily proved using Theorems 3*8 and 3-9 and Remark 3.10. 

1. If X is locally compact, then a motion tt is 

A 

positively Lagrange stable if and only if 
A + (x) is a non-empty compact set. 

2. If a motion tt is positively Lagrange stable, 

X 

then A + (x) is compact and connected. 

3. If a motion tt is positively Lagrange stable, 

then d(xt,A + (x)) 0 as t +°°. 

We are now ready to discuss dynamical systems which are 
noted by their lack of the properties of recursiveness. 

One can easily draw a parallel between these dispersive 
concepts and the recursive concepts just defined. First, 
we need some definitions. 



86 



DEFINITION 4.10 Let x e X. 



1. The motion ir x is said to be positively Lagrange 
unstable whenever cl (y + (x)) is not compact. 

2. The motion is said to be negatively Lagran ge 
unstable whenever cl(y (x)) is not compact. 

3. The motion is said to be Lagrange unstable 

if it is both positively and negatively Lagrange 
unstable . 

4. The point x is called positively Poisson unstable 
whenever x £ A + (x). 

5. The point x is called negatively Poisson unstable 
whenever x £ A~(x). 

6. The point x is called Poisson unstable whenever 
it is both positively and negatively Poisson 
unstable . 

7. The point x is called wandering whenever 
x £ J + (x). 

This is a rather lengthy definition but the definitions 
are what we would expect. We are now ready to define 
dispersive concepts on the whole space X. 

DEFINITION 4.11 The dynamical system (X,R,ir) is said to be 

1. Lagrange unstable if for each x e X, the motion 

it is Lagrange unstable, 

X 

2. Poisson unstable if each x e X is Poisson unstable, 

3. completely unstable if every x e X is wandering. 
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4 . 



dispersive if for every pair of points x,y e X, 
not necessarily distinct, there exist neighbor- 
hoods U x of x and U of y such that U is not 

y x 

positively recursive with respect to U and U 

y y 

is not positively recursive with respect to U . 

X 



One may wonder why such concepts as positively and 
negatively wandering points are not defined. This is 
unnecessary since x e J + (x) if and only if x e J“(x). One 
can see from the definition that a completely unstable 
system must be Poisson unstable since if x £ J + (x), then 
x £ A + (x). The following remark proves the relationship 
between all the sets defined. 



REMARK 4.12 In definition 4.11, each statement implies 
the preceding statement. 



PROOF Suppose (X,R,tt) is dispersive. Then for every pair 
of points x,y e X, there exist neighborhoods U of x and 
U of y such that U is not positively recursive with respect 

y * 

to U . Set x = y. Then there exists a neighborhood of x 
which is not self positively recursive. Hence x is wandering 
so x ?! J + (x) . Since this is true for all x c X, we have 
that X is completely unstable. Next suppose X is completely 
unstable. Then x £ J + (x) for all x e X. Then also 
x £ J"(x) for all x e X. Hence x £ A + (x) and x £ A - (x) 
since A + (x) C J + (x) and A~(x) C J~(x). Thus X is Poisson 
unstable. Finally assume X is Poisson unstable. Then 
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x £ A + (x) and x £ A”(x) for all x e X. Suppose cl (y + (x)) 
is compact. Then by Theorem 3*8, A + (x) is non-empty and 
compact, hence complete. Then by Theorem 4.6 every point 
in A + (x) is non-wandering and by Theorem 4.8, the set P of 
Poisson stable points is dense in A + (x). Thus P is non- 
empty. But if q is a Poisson stable point in A + (x), then 
q e A + (q) and the space is not Poisson unstable. Hence 
we have a contradiction so that cl (y + (x)) is not compact. 

A similar argument shows that cl (y - (x)) is not compact. 
Hence X is Lagrange unstable. // 

The following diagram provides a summary of the ideas 
in the preceding remarks. 




Thus we know that a dispersive system is completely 
unstable, a completely unstable system is Poisson unstable 
and a Poisson unstable system is Lagrange unstable. One 
may wonder if the reverse implications hold. In general 
they do not as the following examples illustrate. 

To show that a Lagrange unstable system may not be 
Poisson unstable we want x e A + (x) but cl (y + (x)) not 
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compact. Since x e A + (x), we know that y + (x) C A + (x) and 
hence cl (y (x)) = A + (x). Thus we must have an example 
where A + (x) is not compact and yet x e A + (x). To this end, 
consider Example 4.3 of the dynamical system defined on the 
torus. We delete the point p from the torus so that the 
torus is no longer compact. Then for any x in the torus 
such that x fi Yg, we have A + (x) = the torus and hence 

+ + J. 

x e A (x) but A (x) is not compact. Thus cl (y (x)) 
is not compact. Clearly for x £ y^ A“(x) = 0 and hence 

cl (y”(x) ) is not compact. Similar results hold for x e y^ 

+ 

and cl (y (x)) is not compact. Thus for all x in the torus, 
cl (y + (x)) is not compact and cl (Y~(x)) is not compact. 

Hence this system is Lagrange unstable but it is not Poisson 
unstable because for any x / y£j x e A + (x). 

To see that a Poisson unstable system may not be 
completely unstable, consider the following phase portrait 
in the x,y-plane. 




The unit circle contains a rest point p and a trajectory y 
such that for each point q e y we have A + (q) = A (q) = (p). 
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All trajectories in the interior of the unit circle have 
the same positive and negative limit sets as y. All tra- 
jectories in the exterior of the unit circle spiral to 
the unit circle as t -*• +«>, so that for each point q in 
the exterior of the unit circle we have A (q) = {p} \J y, 
and A (q) = 0. If we consider the dynamical system obtained 
from this one by deleting the rest point p then the resul- 
tant system is Poisson unstable. However it is not 
completely unstable because for each q e y, J + (q) = y 
so q e J + (q) . 

Finally to see that a completely unstable system may 
not be dispersive consider the following dynamical system 
defined by the differential equations 

x = sin y 
2 

y = cos y 

This system contains trajectories y^ given by 

y k = { (x,y ) : y = kTT + tt/ 2}, k = 0, ±1, ±2,... 

These are lines parallel to the x— axis. Between any two 
consecutive y k 's the trajectories are given by 
y = {(x,y): x + c = sec y}, where c is a constant depending 
on the trajectory. The phase portrait between the lines 
y = _tt/ 2 and y = tt/ 2 is shown in the diagram below. 
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It can be seen that this system is completely unstable 

since for any p e y , , J + (p) = y and for all other points 

-i o 

p not belonging to any y^, J + (p) = 0. Thus p i J(p) for 
any p and hence every point is wandering. However if 
p e y ^ and q e y , then every neighborhood of q is recur- 
sive with respect to any neighborhood of p and the system 
is not dispersive. 

The following result shows that to prove Lagrange 
instability, it is sufficient to prove that the space 
is positively Lagrange unstable. 

THEOREM 4.13 A dynamical system (X,R,tt) is_ Lagrange unstable 
if and only if it_ is_ positively Lagrange unstable . 

PROOF The necessity is trivial. On the other hand, let 
the space be positively Lagrange unstable. Then for each 
x e X, cl (y + (x) ) is not compact. Suppose cl (y~(x)) is 
compact. Then A (x) is non-empty and compact. Since A (x) 
is invariant, yt e A (x) for all t e R and y e A (x). Then 
since A - (x) is closed, cl (y + (y)) C. A (x) and hence 
cl (y + (y ) ) is compact. This is a contradiction. Hence 
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cl ( y" ( x ) ) is not compact for each x e X. Thus the space 
is also negatively Lagrange unstable, yielding the result. // 



For a dynamical system defined in the x,y-plane, the 
following theorem provides an important fact. 

THEOREM *4.1^ For dynamical systems defined by a system of 
differential equations in the euclidean plane , the concepts 
of Lagrange instability and complete instability are 
equivalent . 

PROOF Let this system be given by the differential 
equations 



x = f (x,y) , 
y = g(x,y) , 

where f(x,y) and g(x,y) are continuous functions. 

Clearly if this system is completely unstable, then 

it is Lagrange unstable. So assume the system is Lagrange 

unstable and that there is a point p which is non-wandering. 

Then p e J + (p). The point p cannot be a critical point for 

then cl (y + (p)) would be compact contrary to hypothesis. 

Likewise, there exists a neighborhood V of non-critical 

points about the point p. 

Pi 

Let p = ( ). Since p is not a critical point, then 

p 2 

either f(p 1 ,P 2 ) / 0 or gCp-^P^ ¥ 0. Assume f(p 1 »P 2 ) 5* 0. 
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Since f(x,y) is continuous, there is a neighborhood W about 

p where f(x,y) / 0 for all points ( x ) e V/. Now consider 

*/ 

the function H:R 2 -> R defined by 

H(x,y) = tan ^ ^ f{ " x * ' y ) ^ • Clearly H(x,y) is continuous in 
W. Since H(x,y) = tan ^ [ 7 -] = tan ^(^)> we may consider 

X 

the function H to be a "slope" function since it provides 
us with the slope of any trajectory at a given point. 

Now by the continuity of H(x,y) there is a neighborhood 
U' such that for all points ( X ) in U', 

y 

|H(x,y) = H(p^,p 2 )| < e, where e > 0 is chosen such that 
U' C W. The relationship between the sets V, W, and U' 
is shown in the following diagram. 




Now since |H(x,y) - H(p^,p 2 )| < e for all ( X ) in U', 
we see that the neighborhood U' contains only points whose 
trajectories have slopes that differ from the slope of 
y(p) at p by less than e. Now let 6 > 0 be so chosen 
such that 6 < e and 6 < 1/10. Now let 

U ='{( X ) e U T : | H( x,y ) - H(p 1 ,p 2 )| < <S } . The number 1/10 
was chosen so that we are assured that we can construct 
the following "transversal" through all the trajectories 
lying within U. 
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Let L be a line segment through the point p making 
an angle of 90° with the motion through p. Such a line 
can be constructed since we are working in the plane. Now 
every trajectory in U must cross this line L by our above 
construction of U. 




We claim that each trajectory crossing L in U must do 
so in the same direction as y(p). To see this, let r e U 
and suppose y(r) crosses L in the opposite direction of 
y(p). Let p* be the point where y(r) meets the line L. 
Now let U-L = VJ Ug. It is clear that and are 
disjoint sets. 



Since y(r) and y(p) cross L in different directions, 
we may assume that for some small time t' > 0, p’t’ e 
and pt ’ e U 2 ' Now let L’ be the segment of L which joins 
the points p and p’ and consider irCL'jt’)* Since L’ is a 
connected set, Tr(L',t') is a connected set. Also part 
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of L't' will be in and part of L't' will be in U 2 
since p't' e and pt ' e U 2 . 




Now L't’ is a connected set so there is at least one 
point q such that q e L'H L't'. Thus it follows that 
q = qt ' since it would be impossible that any other 
trajectory passing through L' at time t = 0 would cross 
the point q in time t = t’. For this to happen, we would 
have the following situation. 




Clearly this is an impossibility in U and hence we 
have q =■ qt ' . But this implies that q is a periodic point 
and hence cl (Y + (q)) is compact and we have a contradiction. 
Hence all trajectories cross L in the same direction. 

Now p e J + (p) so there exist sequences {x^} in X 
and (t n ) in R such that x^ -* p and t^ +°° with x^t^ p. 

If any of the x^’s are periodic, then cl (y+Cx^)) is compact 
and we have a contradiction. Hence no x n is periodic and, 
without loss of generality, we may assume each x n lies on 
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L. Now for n sufficiently large, we may assume that 
e ^ since x^t^ -*■ p. Then it follows that there 
exists a T e R such that x fi t (T) lies on L. Call this 
point z. 

It is clear that z lies between the points x^ and p, 
for if not we have the situation shown below and it is 
clear that in this case, x^t^ p is an impossibility. 




Now, let be the segment of L which joins the points 

x^ and z. Since z = x (t + T), we consider the set 
n n n 

(x n [0,t n + T] U L^} = A. This is a Jordan curve and hence 

the interior of A together with A is a compact subset of 
2 

R . Let B denote this compact set. Now p e B so we 
consider y + (p). We claim that pt e B for all t e R + . If 
this is not the case, then y + (p) must cross the set A 
at some point since y + (p) is a connected set. But y + (p) 
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cannot cross A anywhere along x [O.t + T] because different 

n n 

trajectories cannot intersect. Then the only alternative 
is that y + (p) must cross L^. 




But this too is impossible since to do so would mean 
crossing L in the opposite direction from y(x n ). Hence 
pt e B for all t e R + . Hence y + (p) C B and since B is 
compact, we have cl (y + (p))C B. Therefore cl (y + (p)) 
is compact. Thus we have our final contradiction and the 
assumption that p e J + (p) is untenable. Thus the system 
is completely unstable. / 

(The above proof was communicated to me through a private 
conversation with Dr. William Chewning. ) 

We now continue with several characterizations of 
dispersive systems. 

t 

THEOREM 4.15 A dynamical system (X,R,tt) is. dispersive if 
and only if J + (x) = 0 for each x e X. 

PROOF Let the system be dispersive. Assume there is a 
point x e X such that J + (x) 0 ’ 0 • Then if y e J (x) 
there are sequences (x^} anc ^ ^ such that **■ x > 
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+0 ° an< -* X n^n This shows that for any neighborhoods 

U of x and U of y, U t f\ U / 0 as the element x t 
A y x n y n n 

is contained in this intersection. But this contradicts 
the assumption that the system is dispersive because U is 

y 

then positively recursive with respect to U . Hence 

X 

J + (x) = 0 for all x e X. Conversely, let J + (x) = 0 for 
all x e X. We claim that there exist neighborhoods U of 

X 

x and U of y and a T > 0 such that U tfl U =0 for all 
v x y 

t > T. For if not there exist sequences {x } {y } {t } 

— n J n n 

such that x n -*■ x > y n = x n t n , y n -*■ y, and t +», so that 
y t J + (x). But this is a contradiction. Similarly U 

y 

is not positively recursive with respect to U . Hence the 
system is dispersive. // 

THEOREM 4.16 For a given (X,R,ir) the following are equivalent . 
1. The space (X,R,tt) is dispersive . 

^ • For any two points x , y iri X there are neighbor- 
hoods U of x and U of y and a constant T > 0 

x y 

such that U x H U t = 0 for all t , | t | > T. 

3. For any two distinct points x,y in X, y l J + (x) 
and x t J + (y ) . 

PROOF Let (X,R,tt) be dispersive. Then for any two points 

x,y in X, there exist neighborhoods U and U such that U 

x y x 

is not positively recursive with respect to U and U is 

y j 

not positively recursive with respect to U . Hence for some 

T > 0, U x f\ U t = 0 for all t > T and U x t fl U y = 0 for all 

t > T which implies U fi U (-t) = 0. Thus U 0 U t = 0 for 

x y x y 

all t, | t | > T. 
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Assume statement 2 and suppose x e J + (y). Then there 

exist sequences {y^} and {t^} such that y n y and t n •+ +°° 

and y t -*■ x. For n sufficiently large we may assume that 
n n 

v e U and y t e U and t > T. Then we have 
J n y n n x n 

U x n U y t n 0 0 which is a contradiction. The proof that 
y 0 J + (x) is similar. 

Finally assume that for any two distinct points x,y in 
X, y 0 J + (x) and x 0 J + (y). Suppose for all neighborhoods 
U x and U y and any T > 0 we have U x 0 U y t 0 0 for some 
t > T. Since this is true for all neighborhoods, it will 
be true for a nested sequence of neighborhoods about the 
points x and y. Hence consider for each n, the neighborhoods 
S(x,l/n) and S(y,l/n). Now set T = T^ Then 
S(x,l) C\ S(y,l)t 1 0 0 for some t 1 > ' T 1 + 1 - Set T 2 = t l* 

Then there exists a t 2 > T 2 + 2 such that 

S(x, 1/2 ) A S(y,l/2)t 2 0 0 . In this manner we have constructed 
sequences{x n > , (y n ) > (t n > , with x n -»■ x, t n -~*- 00 , ^t^ - y n , 
and y y. V/e note that the process in constructing the 
sequence {t } cannot terminate for if it did, then we have 
found a T. > 0 such that S(x,l/k) A S(y,l/k)t = 0 for all 

rC 

t > T. contrary to our hypothesis. Thus we have x n t n y 

rC 

so y c J + (x) but this is a contradiction and hence U x 
cannot be positively recursive with respect to U y . A 
similar argument holds for U y not positively recursive 
with respect to U . Hence (X,R,tt) is dispersive. // 

A 
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We now introduce the final concept in this section. 

It is that of parallelizable dynamical systems. We are 
headed toward one particular theorem and consequently many 
lemmas and theorems will be presented without proof since 
their proofs are constructive in nature or are not particu- 
larly informative. However, the proofs may be found in 
Bhatia and Szego pages 48 through 53* Formally we begin 
with the definition of a parallelizable dynamical system. 

DEFINITION 4.17 A dynamical system (X,R,tt) is called 
parallelizable if there exists a set S c X and a homeo- 
morphism h: X S*R such that SR = X and h(xt) = (x,t) 
for every x e S and t e R. 

To proceed in the study of parallelizable dynamical 
systems we need to develop a theory of sections. 

DEFINITION 4.18 A set S C X is called a section of (X,R,tt) 
if for each x e X there is a unique time t(x) such that 
xt(x) e S. 

The set S in Definition 4.17 is in fact a section of 
(X,R,tt) for if z is any element in X, then z = xt for some 
x e S and t e R. If this was not the case, then SR ^ X. 

Then the value (-t) is precisely t(z) given in Definition 4.18. 

It is clear that zt(z) = z(-t) = x e S. To see that (-t) 

is unique, assume there exists (-t 1 ) e R such that 

z(-t') = y e S. Then xt = yt*. Now using the homeomorphism 
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h, we have h(xt) = h(yt') which implies (x,t) = (y,t'). 

Thus we have x = y and t = t'. Thus x(z) is unique. 

Then to summarize, it is clear that if a dynamical 
system is parallelizable it must have a subset S which 
contains a point of every trajectory so that SR = X. More- 
over, the set S can contain only one point of every trajec- 
tory for otherwise t(x) would not be unique. To see that 
this last statement is true, suppose there is a trajectory 
which intersects the section twice. Then there are points 
x,y in this trajectory such that x,y e S and xt = y for 
some t > 0. Then t(x) = 0 since xo = x e S but also 
xt = y e S and hence t = x(x). Since t > 0, this contradicts 
the uniqueness of t(x). This implies in particular that 
the dynamical system cannot have any critical points or 
periodic trajectories. To see how a parallelizable system 
compares with the other dispersive concepts defined in this 
section, consider the following example. 

2 

EXAMPLE 4.19 Consider a dynamical system defined in R by 
the differential equations 



dx 

dt 



f (x,y) , 



dy. = 

dt 



0 , 



where f(x,y) is continuous, and moreover f(x,y) - 0 whenever 
the point (*) is of the form with n a positive integer. 
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For simplicity we assume that f(x,y) > 0 for all other 
points. This system is sketched in the diagram below. 












/>L _ 




r vl ** 

>* 




i — ? m — : 




! s** ■ S" ll\ 

. - . yj • — > — 




^ v-> ^ 

> — 



■ 11 1 " /gt ? 1 * 

Now we form the system in which we are interested by 
deleting the sets 



In = ^ ^ y ^ : x - n ’ y = 1//n ^ » n = 1 , 2,3,... 



from the plane R 2 . This system is sketched below. 












”X 



It is clear that this system is dispersive since the 
only way two neighborhoods U x of x and U y of y could be 
positively recursive with respect to each other would be 
if the points x and y lie on the same trajectory. Then 
clearly we could find a T > 0 for which U x H U y t = 0 or 
U A U t = 0 for all jt | > T. Hence this system is disper- 
sive but it is not parallelizable since no section S can 



exist for which t(x) would be continuous. This fact can 
be seen when it is realized that the only section possible 
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is one that looks like y = 1/x + c for any constant c 
together with some section which intersects the trajectories 
below the x-axis. If this "lower" section intersects the 
trajectory on the x-axis, then t(x) could not be continuous. 

If the lower section looks like y = -1/x + c, then t(x) 
would be continuous but no section will intersect the tra- 
jectory on the x— axis . Hence the system is not parallelizable . 

We can get an example of a parallelizable dynamical 
system from the previous example by considering the 
differential equations 

dx _ 
dt * 



This system is sketched below. 





W 














' ( 5 ) 


1 







Clearly any line parallel to the y-axis is a section of 
the system. 

These examples also give an intuitive feeling for the 
notion of a parallelizable dynamical system in that all 
trajectories are in some sense parallel. 
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The function t(x) will be basic in our development. 

In general the function t(x) need not be continuous however 
for a dynamical system to be parallelizable there must 
exist a section where t(x) is continuous on X. This fact 
is proved in the following theorem. 

THEOREM 4.20 A dynamical system (X,R,tt) is_ parallelizable 
if and only if it has a section S with x(x) continuous on X. 

PROOF Suppose (XjRjTt) has a section S with t(x) continuous 
on X. Then clearly SR = X. We now define the homeomorphism 
h: X *»■ SxR by h(x) = (xx(x) ,-x(x) ) . To see that h is one- 
to-one let h(x) = h(y) for x and y in X. Then xt(x) = yx(y) 
and -x(x) = -x(y). Hence xx(x) = yx(x) e S which implies 
that x = y for otherwise x(y) is not unique. Hence h is 
one-to-one. The function h is continuous by the continuity 
of x ( x ) and the continuity axiom. The inverse 
h" 1 : SxR -v X is given by h _1 (x,t) = xt and is clearly 
one— to— one and continuous. Thus h is a homeomorphism so 
the definition of a parallelizable dynamical system is 
satisfied. 

Now assume (X,R } tt) is parallelizable. Then the set S 
in the definition is a section of X. Since for any x e X, 
x = yt for some y e S and t e R, we set x(x) = — t. Then 
xx ( x ) = x ( — t ) = y e S. The continuity of x(x) follows 
from the continuity of h in the definition. // 
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To continue this development we need to define special 
types of sections. This is done in the following definition. 



DEFINITION 4.21 An open set U in X will be called a tube 
if there exists a x > 0 and a subset S C U such that 
SI C U and for each x e U there is a unique x(x), 

T 

| x ( x) | < x, such that xx(x) e S. Here I T = (-x,x). 
is also called a x- tube with section S, and S a (x - U)- 
section of the tube U. If I T = R then U is an °°-tube and 
S an (°° - U)-section. 

DEFINITION 4.22 Given an open “-tube U with a section S 
and x ( x ) continuous on U, and given sets N and K such that 
N C K C S where N is open in S and K is compact, we shall 
call KR the compactly based tube over K. Then indeed x(x) 
restricted to KR is continuous on KR. 

Proceeding toward our final result, the next theorems 
and lemma are presented without proof. 

THEOREM If X is locally compact and separable, and if 
every x e X .is a wandering point , then there gxis. _ s a 
countable covering {K n R} of X by compactly based tubes 
K^R each with x^Cx) continuous on K n R. 

LEMMA A compactly based “-tube U with sectio n K of a 
dispersive dynamical system (X,R,tt) is close d in X. 

We also need the following important topological 
theorem. 
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TIETZE EXTENSION THEOREM Let X be a normal topological 



space , let A be a closed subset of X and let f be. a 
continuous function on A to the closed interval [a,b]. 

Then f has a continuous extension g which carries X into 

[a,b] . 

The final lemma needed to prove our result is quite 
instructive and is completely presented below. 

LEMMA Let , U 2 be_ two compactly based tubes of a dlsper - 
sive dynamical system with sections K^, K 2 and continuous 
functions t 1 (x) and x 2 (x) respectively . If D u 2 ^ 0 
then U = U U 2 is a compactly based tube with a section 
K D and a continuous function x(x). Moreover, 
time distance between and K 2 along trajectories in 
U A U 2 is less than x ( > 0 ) , the time distan ce between K 
and K 2 along tra.j ectories in U is also less, than x . 

PROOF In this proof, the following diagram will be useful. 




By the previous lemma, U 1 and U 2 are closed. They are 
invariant because each tube contains a section K 1 and K 2 
respectively and yc ^ and KjRC U r Hence U-^ A U 2 is 
invariant and closed. Further, K 2 fl U 1 is compact and 
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non-empty, as shown in the diagram. Set S 2 = D 
and f| U 2 . Any trajectory in 0 U 2 intersects 

in exactly one point and intersects S 2 in exactly one 
point. Thus for any x e 0 U 2 , t^(x) = t 2 (x) + t^(xt 2 (x)). 

This is true because xt^Cx) = xt 2 (x) (t^(xt 2 (x)))= x(t 2 (x) + 
Xi(xt 2 (x))) and there are no rest points or periodic 
trajectories in a dispersive system. The function is 
continuous on S 2 which is compact and now we apply the 
Tietze extension theorem. This is possible since every 
metric space is normal. Thus the continuous function x^ 
is extended to a continuous function t defined on K 2 where 
T (x) = t 1 (x) for x e S 2 - If t 1 (x) e for x e S 2 , 

we have t(x) e for x e K 2 - 

Notice now that {xx(x): x e S 2 > = S.^ and t(x) being 
continuous, {xt(x): xe K 2 ) is compact as K 2 is compact. 

Set K = K 1 U (xt (x) : x e K 2 > and define x*(x) on 
KR = K X R U K 2 R as follows: 



T *(x) is continuous on KR and we need only verify that if 





already been proved since t(x) t^(x) for x e S 2 « 
the lemma is proved. 



Thus 



// 
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Finally, we prove the theorem that is the main focus 
of this section on parallelizable dynamical systems. 

THEOREM 4.23 A dynamical system (X,R,tt) on a locally 
compact separable metric space X is parallelizable if and 
only if it is dispersive . 

PROOF Let the dynamical system (X,R,tt) be dispersive. By 

Theorem 4.20 it is sufficient to prove that X has a section 

S with t(x) continuous on X. By a previous theorem, there 

is a countable covering (U n > of X by compactly based tubes 

U with sections K and continuous functions x (x). VJe 
n n n 

replace this covering by a like covering (U n ) of compactly 

1 

based tubes which we construct as follows. Set = K , 

and U 1 = U 1 . Beginning with U 1 and U 2 we use the previous 

1 2 

lemma to enlarge K to a compact set K , thus obtaining 

pi 2 

the compactly based tube U = U U U 2 with x (x) continuous 

o n 

on U . This leaves K 1 unaltered. Having found U in the 

same manner, we take it together with U n+ -^ and construct 

U n+1 with K n+1 Z>K n , and x n+1 (x) continuous on U n+1 . Now 

set S = VJ K n , then X = SR and the function x(x) defined by 

x(x) = x n (x) for x e U n is continuous on X, with the property 

that xx ( x ) e S. Moreover x(x) is unique for each x e X. 

Thus X has a section S with continuous x(x) defined on X. 

The system (X,R,tt) is thus parallelizable. 

Conversely, suppose the system (X,R,tt) is parallelizable. 

Then there exists a section S with x(x) continuous or. X. 
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Suppose J + (x) 0 0 for some x e X. Let y e J + (x). 

Then since y e X, there exists y e S and F e R such that 

y 

yt - y. Now since J + (x) is invariant, y(-F ) e J + (x) 

* y 

which implies that y e J (x). ‘Then there exist sequences 

{x^} and {t n ) with x^ *► x, t n •* +~, and x^t^ y. Nov/ 

for each x , there exists x e S and F e R such that 
n n n 

x F = x . Similarly there exists x e S and F e R such 
n n n x 

that xt = x. Then since x t -*■ y. it follows that 
x n n 

x F (t ) -*■ y. 
n n n 

Consider for each fixed k, k = 1, 2, the sequence 

{x F (t, )}. Since x t -> xt , we have by the continuity 
nn k nn x* 

axiom that x F (t, ) -*■ xF (t, ). Hence we may assume by 
n n k x k 

taking subsequences if necessary that 

d(x F (t, ),xF (t )) <_ 1/k for k = 1, 2, ... Then it follows 
n n k x k 

that xF (t ) -> y since 
x n 

d( y>* F x (t n >) i d( y> 7 n 5 'n (t n )) + d( *nW > - 



d(y,x t (t )) + 1/n. Now since xt (t ) y, let 
n n n x n 

t = t + t . Since t ■+ +°°, it follows that x n +°°. 
n x n n * n 

Hence we have that y e A + (x). 




Now for each n set y n = xT n * Then y n + y • Since 

we may assume that t > n for each n. Let U be an 

n y 



+ CO 



> 
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arbitrarily small neighborhood of y. Then since y f ->■ y } 
y n e for all n > N where N is a large enough integer. 




We know that t(x) is continuous and since i(y) = 0, it 

follows that x(y ) is close to 0 for all n > N. Hence we 

may assume that |x (y )| < 1/n for each n > N. But this is 

impossible since T (y n ) is unique and we have already 

determined that x(y ) =-x < -n for all n. Thus we have 

n n 

a contradiction so it is impossible that y e J + (x) and 
hence J + (x) = 0 for all x e X. Then by Theorem 4.15 we 
have that (XjRjir) is dispersive. // 

The relationships bet ween the various dispersive concepts 
is now conveniently summarized below. 




Ill 



V. STABILITY THEORY 



In this section we will study stability theory first 
from the abstract point of view that we have established 
and then from an applied point of vi ew that uses systems 
of differential equations and studies the eigenvalues of 
such a system. To begin we need some basic definitions. 
Throughout this section , it will be assumed that the space 
X is^ locally compact and the set M is a non - empty compact 
subset of X unless explicitly stated otherwise . 

DEFINITION 5.1 With a given M C X we have 

1. the set A (M) = {x e X: A + (x)f\ M ¥ 0), 

w 

2. the set A(M) = {x e X: A + (x) 0 0 and A + (x) C M} . 

The set A (M) is called the region of weak attraction and 

w — — 

the set A(M) is called the region of attraction of the set 
M. Moreover, any point in A W (M) or in A(M) is said to be 
weakly attracted or attracted to M respectively. 

From this definition it is clear that a point x is 

weakly attracted to a set M if and only if there is a 

sequence {t } in R with t -*■ +°° and d(xt ,M) -*■ 0. Also a 
n n n 

point x is attracted to a set M if and only if d(xt,M) 0 
as t + +°°. 

The above attractor sets are important in the development 
of stability theory and the following theorem provides 
needed information about them. The following lemma is useful. 
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LEMMA 5.2 Let X be_ any metric space and let x e X. Then 
for every t e R, 

1. the_ set A + (x) = A + (x)t = A + (xt), 

2. the set_ J + (x) = J + (x)t = J + (xt). 

PROOF The first equality in statements 1 and 2 follows 
from the invariance of the sets A*(x) and J + (x) respectively. 
To see the second equalities, consider the case 

•f> i . 

A ( x ) t = A (xt). Let z e A (x)t. Then there is a y e A (x) 

with z = yt and a sequence (t }, t -*• +°° with xt -*• y. 

n ’ n 3 n 

Then by the continuity axiom, xt n (t) -> yt. But 

xt (t) = xt(t ) and since t ->+<», we must have that 
n n n 

+ + + 

yt e A (xt). Thus A (x)t <Z A (xt). The argument is clearly 
reversible so we also have that A (xt)C A (x)t. Hence 
statement 1 is proved. The proof of the second inequality 
in statement 2 is entirely analogous. // 

THEOREM 5* 3 For any given M, A VJ (M) ZD A(M) and the sets 
A., (M) and A(M) are invariant. 

W r 1 “ " 1 

PROOF Let x e A(M). Then A + (x) f 0 and A + (x) C M. Then 

clearly A + (x) 0 M ¥ 0. Hence x e A (M) . To see the 

w 

invariance of the sets A (M) and A(M) , consider the case 

w 

for A (M). Let x e A,,(M) and let t e R. Now 
w w 

4* 4. 4- 

A (x) H M j- 0 and from Lemma 5-2, A (x) = A (xt). Hence 
A + (xt) n M / 0 which implies that xt e A (M). Since t e R 
was arbitrary, we have that A (M) is invariant. The proof 
of the invariance of A(M) is entirely analogous. // 
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Using these sets, we are now ready to define the various 
types of stability to be investigated in this section. 

DEFINITION 5.4 A given set M is said to be 



1. a weak attractor if A„(M) is a neighborhood of M, 

2. an attractor if A(M) is a neighborhood of M, 

3. stable is every neighborhood U of M contains a 
positively invariant neighborhood V of M, 

4. asymptoti cally stable if it is stable and is 
an attractor, 

5. unstable if it is not stable. 



Geometrically, M is stable if given any neighborhood U 
of M, there is a neighborhood VC U such that any trajectory 
that is within V at time t = 0 remains within V for all 
t > 0. If, in addition, all trajectories approach M as 
t increases, then M is asymptotically stable. Finally 
M is unstable if every neighborhood of M contains trajec- 
tories which are arbitrarily close to M at time t = 0 but 
which eventually leave the neighborhood for all time. These 
three cases are illustrated in the diagram below. 
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We have seen specific instances of these cases in 
several of our examples. For instance in the pendulum 
problem, (Example 2.1), the point (^) is stable; in the 
damped spring problem (Example 2.2), the point (q) is 
asymptotically stable and in Example 3*2, if the system is 
restricted to the interior of the unit circle, the point 
(q) is unstable. In this latter example, if x is any point 
on the unit circle, then the point x is a weak attractor 
although it is not stable. In order to see an example of 
an attractor which is not stable, consider the following 
example : 

Let the planar dynamical system be defined by the 
differential equations (in polar coordinates): 

r = r(l - r) , 

9 = sin 2 (6/2). 

\ 

The trajectories of this system are sketched below. They 
consist of two critical points at (q) and (q), & trajectory 
on the unit circle with (g) as the positive and negative 
limit set of all points on the unit circle. All other points 
have the point (q) as their positive limit set and hence 
the point (J) is an attractor. It is not stable because 
every neighborhood will contain part of the unit circle 
where 0 > 0. Hence no positively invariant subset of this 
neighborhood can exist. 
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We are now ready to prove some theorems about the 
concepts just introduced. 

THEOREM 5.5 If M is a weak attractor ( respectively , an 
attractor ) , then the set A w (M) ( respectively , A(M)) is_ an 
open neighborhood of M. 

PROOF Let N denote either of the sets A W (M) or A(M) . Then 
N is an invariant neighborhood of M. Consequently, bdy N 
is invariant by Theorem 2.10, is disjoint from M, and is 
closed. Then for each x e N, A + (x) 0 M / 0, whereas for 
each x e bdy N, A + (x) C bdy N since bdy N is invariant 
and closed. Since bdy N A M = 0, we conclude that 
N 0 bdy N = 0. Thus N is open. // 

We note that if a set M is stable then it is the inter- 
section of positively invariant neighborhoods and M is 
therefore positively invariant. Immediately from this we 
see that if the singleton (x) is_ stable , then (x } must be 
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positively invariant so consequently , x must be a critical 
point . 

The next two theorems show different characterizations 
of stable sets but first we establish the following lemma 
and its corollary which do not depend on the local compact- 
ness of X. 

LEMMA 5 • 6 Let X be an arbitrary metric space . Let x e X 

+ + + 

and w e A ( x ) . Then J ( x ) C J ( w ) . 

PROOF Given w e A (x) and any y e J (x), there exist 
sequences (xM , -*■ +°°, xx^+w, and {t^} and (x n ), 

x -*■ x. t ' +°°, x t 1 y. We may assume, if necessary by 

choosing subsequences, that t^ - x^ > n for each n. 

Consider for each fixed k, k = 1, 2, ..., the sequence 
{x n Tjp- By the continuity axiom x n T k"*" XT k> k = 1 5 2 > 

We may, therefore, without loss of generality, assume that 
for each fixed k, dtxx^x^) 1 1A for n > k. ' This shows 
that xt’ w, because d(w,x x')<d(w,xx’) +d(xx' ,x t') <. 
d(w,xx^) + 1/n. Now notice that x^ = x^Ct^ -x^), 

and x n t I ! 1 y, x n T 1 l 1 * w > anc * ~ T n > n ’ Hence 

y e J + (w). As y e J + (x) was arbitrary, we have 

J + (x) C J + (w), and the lemma is proved. // 

COROLLARY Given M and x e A W (M) , then J (x) J (M) C D (M) . 

PROOF Indeed for any x e A W (M), A + (x) 0 M / 0. Let w 
be any element in A + (x) C\ M. Then by the previous lemma, 
J + (x) C J + (w). But J + (w)C J + (M) and we know that 
J + (M) C D + (M). Hence J + (x)C J + (M) C. D + (M) . // 
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Vie now give an important characterization of stability 
of a set M. 

THEOREM 5.7 A set M is stable if and only if D + (M) = M. 

PROOF Let D + (M) = M and suppose if possible that M is not 
stable. Then there is an e > 0, a sequence {x^}, and a 
sequence {t n }, with t fi > 0, d(x n ,M) -> 0, and d(x t ,M) > e . 

We may assume without loss of generality that e > 0 has 
been chosen so small that S[M,e] and hence H(M,e) is compact 
where H(M,e) = {x e X: d(x,M) = e}. This is possible since 
X is locally compact. Further, we may assume that 
x x e M. Vie can now choose a sequence {t},0<t <t, 

such that x^x n e H(M,e), n = 1, 2, .... Since H(M,e) is 
compact, we may assume that x t y e H(M,e). Then clearly 
y e D + (x)d D + (M), but y t M. This contradiction shows that 
M is stable. Conversely, assume that M is stable. Then 
given any neighborhood U of M there is a positively invariant 

neighborhood V of M with VC U. Since for any x e M, 

+ + 

D (x) C cl WR for any neighborhood VI of x, we get 
D + (x) C cl V since V is positively invariant. Thus 
D + (M) C cl U for any neighborhood U of M. Hence 
D + (M)C n(cl U: U is a neighborhood of M} = M as M is 
compact. Since M C. D + (M) always holds, we have D + (M) = M. // 

The following remark shows what results when we combine 
several of the concepts under discussion. 
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REMARK If M is_ stable and Is a weak attractor , then M is_ 
an attractor and consequently asymptotically stable . 

PROOF Since M Is a weak attractor, A w (M) is a neighborhood 
of M. We need to show that A w (M)C A(M). Let x e A v/ (M) . 

Then A + (x) 0 0 since A + (x) H M 0 0 and A + (x) C. J + (x) <L D + (M) 
by the previous corollary. Since M is stable we have 
D + (M) = M by Theorem 5*7 and consequently A + (x) C M. Thus 
x e A(M) and A(M) is a neighborhood of M. Thus M is an 
attractor. Since M is stable and is an attractor, M is 
asymptotically stable. 

The various implications that have been proved thus far 
are summarized in the following diagram. 




We would now like to confine our attention to the study 
of stability for systems of differential equations, he 
shall begin by considering autonomous systems of constant 
coefficient linear equations; that is, systems of the form 



X 



AX 
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where A is an n*n real valued matrix. The origin in R n 
is a critical point for such a system, and when zero is 
not an eigenvalue for A (so that is when A is non-singular) 
the origin is the only critical point. The solutions of 
this system can be described in terms of the eigenvalues 
of A in the following manner. 

Let X = a + bi be a complex eigenvalue for the nxn 
real matrix A, and let be an eigenvector in C n belonging 
to X. Then the functions 



X^(t) = e at (G^cos bt + H^sin bt), 
X 2 (t) = e at (H x cos bt - G^sin bt), 



are linearly independent solutions of X = AX, where 



G 



X 



E x + E x 

2 




^ E X " V 
2 



and 



[Kreider, Kuller, Ostberg, 

p. 255] 



Now when a ± bi, b > 0, are complex conjugate eigenvalues 
of multiplicity m, this system has 2m linearly independent 
solutions constructed from functions of the form 



n at . k 
G x e t 



cos bt 



and 
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„ at, k . , , 

H^e t sin bt, 

where k is an integer such that 0 £ k £ m - 1. Then 
assuming t > 0, 



| |G x eattk c ° s bt | | <_ | [ | | e at t k , 

and since | | | | is a positive constant, 

lim | | G e at t k cos bt | | =0 if a < 0. 
t-*-<» A 

On the other hand ||G x e at t k cos bt | | is unbounded if a > 0. 
Similarly , 

lim | |ihe a ^t k sin bt | | = 0 if a < 0, 
t-**> A 

but ||H^e a ^t k sin bt | | is unbounded if a > 0. 

Thus every traj ectory of this system arising from a 
pair of complex conj ugate eigenvalues for A will approach 
the origin of R n as_ t +» if the real part of these eigen - 
values is negative ; it will depart arbitrarily far from 
the origin , if the real part of these eigenvalues is positive . 
In the first case, the trajectories exhibit the properties 
required for asymptotic stability whereas in the second 
case they exhibit instability. This reasoning clearly 
applies equally well to trajectories arising from real 
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eigenvalues for A. It remains to consider the case where 
A admits a pair of pure imaginary eigenvalues ± bi. If 
these eigenvalues have multiplicity one, the corresponding 
solutions of the system are _ constructed from the functions 

G^cos bt and H^sin bt 

and the fact that cos bt and sin bt are bounded as 
t +co yet do not tend to zero implies that the resulting 
trajectories have the properties required for stability, 
but not asymptotic stability. On the other hand, if the 
multiplicity of ±bi is greater than one, the system has 
solutions involving functions of the form 

G^t^cos bt and H^t^sin bt 

with k >_ 1. When this happens, the origin is unstable since 
Jt^cos bt | and |t^sin bt | 

are unbounded. Thus we have proved the following theorem. 

* 

THEOREM 5*8 If X = AX is_ an nxn linear autonomous system 
whose coefficient matrix is nonsingular , then the origin 
in R n is_ 

1. asymptotically stable if the real parts of all 
of the eigenvalues of A are negative ; 
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2. 



stable , but not asymptotically stable , 1 f A 
has at least one pair of pure Imaginary eigen - 
values of multiplicity one , no pure imaginary 
eigenvalues of multiplicity exceeding one, and 
no eigenvalues with positive real parts ; 

3. unstable otherwise . 

Consider the following plane autonomous systems: 

• • • 

1. x = y 2 . x = y 3- x = -x 

• • • 

y = x y = -x y=-y 

The coefficient matrix for the first system is 



A = 



r- 

0 

1 



1 

0 



with eigenvalues X = ±1. Hence the origin is unstable 
since one of the eigenvalues has a positive real part. 
For the second system, 




1 

0 



This system is precisely that of Example 2.3 and has 
eigenvalues X = ±i. Thus the origin is stable but not 
asymptotically stable. This is exactly as we would expect 
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by recalling the geometry of Example 2.3. 
For the third system, 



-1 



A = 



0 



0 




This system has X = -1 as its only eigenvalue and thus 
the origin is asymptotically stable. 

The above results can be used to obtain a complete 
description of the trajectories of a linear plane autonomous 
system of the form 



x = ax + by 
y = cx + dy 



where the matrix 




«n 



is nonsingular. This description depends upon the nature 
of the eigenvalues of A which are the roots of the 

characteristic equation 



X 2 - (a + d)X + (ad - be) = 0. 
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Since the discriminant of this equation is given by 



A " (a + d) 2 - *1 (ad - be) = (a - d) 2 + fJbc, 

it follows that A^ and X 2 will be 

1. real and distinct if A > 0, 

2. real and equal if A = 0, 

3. complex conjugates if A < 0. 

Moreover, when A^ and A 2 are the roots of this system, we 
have (A - A^)(A - A 2 ) = 0. Thus we have 

A 2 - (A 1 + A 2 )A + A x A 2 = 0, 



so 



A^ + A 2 = a + d, A-^A 2 = a< ^ - be = det A 

Thus when A >_ 0, A-^ and A 2 will have the same sign if and 
only if ad - be > 0, and will then be positive or negative 
according as a + d is positive or negative. On the other 
hand, if A < 0, A^ = a + pi, A 2 = a - pi, P > 0, and 
a + d = 2a. 

In this case A 1 and A 2 will be pure imaginary if and only 
if a + d = 0. Otherwise they will have a nonzero real part 
which agrees in sign with a + d. 

Now, combining these observations with Theorem 5-8, we 
have the following description of the origin for a plane 
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autonomous linear system with constant coefficients. 



1 . 



2 . 



3 . 



A 



A 



A 



< 



> 



0: 



0 : 



0 : 



! 

! 



a + d < 
a + d > 
a + d = 
ad - be 
ad - be 
ad - be 
ad - be 
ad - be 



0 

0 

0 

> 

> 

< 

> 

< 



0 

0 

0 

0 

0 



asymptotically stable 

unstable 

stable 

and a + d < 0 asymptotically 

stable 

and a + d > 0 unstable 

unstable 

and a + d < 0 asymptotically 

stable 

unstable 



Thus far we have not discussed the cases where A is 
singular. These cases are rather uninteresting but for 
completeness, they are presented here. 

Suppose A is the 2x2 zero matrix. Then X = 0 for all 

2 2 
X e R which implies that each point in R is a critical 

point and consequently each point is stable since every 
point within some neighborhood of a point, remains within 
that neighborhood as t •+ +«°. 

Suppose A is a 2x2 matrix with zero as its only eigen- 
value. Then A must be of the form 



A = 



a 

a 



-a 



-a 



where a is some real number. 
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To find the eigenvectors X Q associated with A, we set 
AXq = 0. Then all eigenvectors are critical points and 
are given by the equation 




which implies that x^ = x 2 * Thus all points on the line 
through the origin with slope = 1, are critical points. 
Now 

X 1 ~ aX l ~ aX 2 
x 2 = ax 1 - ax 2 



so we have x^ = x 2 which implies that x 2 = x 1 + c where c 
is any constant. Thus the remaining trajectories also have 
slope = 1, and intersect the x 2 ~axis depending on c. 

Their directions are shown in the sketch below. It is clear 
that all critical points are unstable. 
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Thus far we have studied only linear systems. To 
study non-linear systems we need the following idea of a 
linear approximation. Consider the system 

X = F(X) 

where X is a vector in R n . We shall show that whenever F 
is a function of class (that is, F has continuous first 
partial derivatives) in a region ft of R n containing the 
origin, then we may replace F by a function of the form 

JX + G(X) , 

where J is an nxn matrix and G(X) is "small" in comparison 
with JX when ||x|| is small. Precisely, we state 

DEFINITION 5*9 Let F be as above and let Xq be a point in 
ft. Then a function L = L(X) is said to be a linear 
approximation to F at Xq if L is linear on R n and 

CF(X) - F(X 0 )] - [L(X) - L(Xq)] 

lim — = 0. 

I |X-X 0 | 1-0 1 1 X - X 0 1 1 

Moreover, the existence and uniqueness of the linear 
approximation can be guaranteed if F is of class c\ 
[Ostberg, Kreider, Kuller pp . *J17, ^ 1 8 ] . This linear 
approximation is in fact the n*n Jacobian matrix of F at 



128 



every point X Q . Specifically, 



J = 



3x n 



3F 



n 



dX- 



3x 



n 



3F 



n 



3x 



n 



X = X, 



where 



F(X) = 



F 1 (X) 



F (X) 
n 



and F ± :R n R for all i. 

To gain a better understanding of the linear approxima- 
tion, we give some examples. 



We find the linear approximation about (q) for each of 
the following systems: 

1. x = P^Xjy) 

y = P 2 (x,y) 



where and P 2 are polynomials with P^(0,0) = P 2 (0,0) = 0. 
Then F^(x,y) = P^(x,y) and F 2 (x,y) = P 2 (x,y). Now if 



129 



P-^(x,y) a^x + b^y + c^xy + a2X^ + b2y^ + higher-order terms, 
P 2 (x >y ) = d^x + e^y + f^xy + d2X^ + e ^ y^ + higher-order terms, 



we have 



= a ]_ + c 1 y + 2a2X + higher-order terms and 

8F 1 

= b.^ + c-^x + 2 bgy + higher-order terms 
Now evaluating these partials at the point (q) , we have 

3F-, 3F 

— ‘ a i and sr ’ b i 

Similar results hold for the partials of F^. Thus we 
replace F(X) by JX + G(X) where 









- 2 . 2 , n 




a i b i 




c 1 xy + a 2 x + b 2 y + • • • 


J = 




, and G(X) = 






d i e i 

«■ 




2 2 

f^xy + d 2 x + e 2 y + • • • 



Hence the linear approximation is given by the following 
system: 

X = JX. 
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2. x = sin x + e y - 1 
y = xy 

Now, sin x and e y have the following Taylor series 
expansions : 



sin x 




e 



y 



2 

i + y + |r + • • • • 



Hence we have 




x + 



y " 3T 



2 

y 

— + higher-order terms 



y 



xy 



and applying the procedure used in the previous example, we 
have that the approximating system is given by 



X = 



JX 



where 




1 

0 

J 



To analyze this system as to stability, we note that the 
matrix J is singular and hence our theorem does not apply. 
Thus we must analyze J using the same method we used in 
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the two previous cases of singular matrices. We see that 
J has eigenvalues X = 0,1 and hence does not fit into the 
category of the other two examples. However, it is an 
easy task to determine that all eigenvectors lie on the 
line of slope -1 through the origin and are all critical 
points. Prom the approximating system, it is easily seen 
that y = 0 and hence all trajectories not on the line of 
critical points are parallel to the x-axis. Their direc- 
tions are easily determined from the equations. They are 
shown in the following diagram. 




It is clear that each critical point is unstable. 

One might question now whether the complete theory we 
developed for linear systems also applies to the approxima- 
ting systems for non-linear systems. Such is the case for 
asymptotically stable and unstable systems. That is, i- 
the linear approximating system is unstable, then the non- 
linear system is unstable. However, no definite conclusion 
can be made about stable systems. It is recalled that for 
a system to be stable, it could have no eigenvalues with 
positive real parts. A non-linear system may have 
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eigenvalues with small positive real parts which do not 
appear in its approximating system. These facts are 
proved in the next section and are presented here Just 
for information. 
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VI. LIAPUNOV STABILITY 



The basic feature of stability theory a la Liapunov is 
that one seeks to characterize stability and asymptotic 
stability of a given set in terms of a non-negative real- 
valued function defined on a neighborhood of the given set. 
This technique was devised by the Russian mathematician 
A. A. Liapunov, and is known as Liapunov's direct or second 
method. It is based on the well-known fact that a physical 
system loses potential energy in a neighborhood of a point 
of stable equilibrium. More precisely, a point of stable 
equilibrium for a physical system is a point at which the 
potential energy of the system has a local minimum. This 
fact is known in physics as Lagrange's theorem. We proceed 
with the following lemma. 

LEMMA 6 . 1 Let the phase space X be arbitrary and let K C X. 
Let <J> be_ any continuous real - valued function defined on K 
such that cf>(xt) £ $(x) whenever x[0,t] C K, t £ 0 . Then 
if for some x, cl (y + (x)) C K, we_ have 4>(y) = 4>(z) for 
every y,z e A + (x). 

PROOF Now there are sequences {t n } and {x^} in R such 

that t -^ +<», x +oo, and xt -* y, xx -*■ z. We may assume 

by taking a subsequence that > t^ for each n. Then 

clearly <|>(xt ) £ cKxt^) since xx^ = xt n ( T n “ t n )» 

t - t >0, and x t [0,t - t ] <^ K. Thus proceeding to 

n n * n n n n 



134 



the limit we have by the continuity of 4>, <})(y) >_ <p(z) . 

A similar argument shows that 4> ( 2 ) > 4> ( y ) . Thus the lemma 
is proved. // 

The best known result on asymptotic stability is the 
following. 

THEOREM 6.2 A compact set M C X is_ asymptotically stable 
if and only if there exists a continuous real - valued 
function E defined on a neighborhood N o£ M such that 

1 . if x e M , we_ have E ( x ) = 0 and if x l M , we 
have E(x) > 0, 

2. for x l M, t > 0 and x[0,t] C N, we have 
E(xt ) < E(x) . 

PROOF Assume that a function E as required is given. 

Choose a > 0 such that S[M,a] C N and is compact. Let 
m = min{E(x):x e H(M, a) } where H(M,a) = {x e X:_d(M,x) = a}. 
Such a minimum exists because E is continuous and H(M,a) 
is a closed and bounded set. By statement 1 and the 
continuity of E we have m > 0. Set 

K = {x e S[M,a]:E(x) < m>. Then K is compact and because 
of statement 2, K is positively invariant. This establishes 
that M is stable as K is a positively invariant neighborhood 
of M. To see that M is an attractor, choose any compact 
positively invariant neighborhood K of M with KC N . Then 
for any x e K, y + (x) C K and since K is compact, 
cl (y + (x))C K. Hence cl (y + (x)) is compact and we have 
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0 ¥ A + (x) C K, and Lemma 6.1 shows that E is constant on 
A (x). But this shows by statement 2 that A + (x)C M. Thus 
M is an attractor and, consequently, asymptotically stable. 

Conversely, let M be asymptotically stable and A ( M ) 
its region of attraction. For each x c A(M) define 

4 >(x) = sup{d(xt,M): t >_ 0}. 

Indeed <f>(x) is defined for each x c A(M) because if 
d(x,M) = a, then there is a T > 0 with x[T,+oo) C. S(M,a). 

This is because x e A(M). Thus 

<p(x) = sup{d(xt,M): 0 <_ t <_ T} . 

Since d(xt,M) is a continuous function of t, 4>(x) is defined. 
Now ( x ) has the following properties: <f>(x) = 0 for x e M, 
<f>(x) > 0 for x i M, and ( xt ) <_ <J>(x) for t >_ 0. This is 
clear when we remember that M is stable hence positively 
invariant, and that A(M) is invariant. Thus if <f>(x) is 
defined for any x e A(M) it is defined for all xt with 
t e R + . We further claim that this ( x ) is continuous in 
A(M). Stability of M implies continuity of 4>(x) on M. 

For x 0 M, this proof requires a knowledge of uniform 
attractors and it will not be presented here. However, 
this function is continuous in A(M) [Bhatia, Szego p. 67]. 
This function may not be strictly decreasing along 
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trajectories in A (M) which are not in M and so rnay not 
satisfy statement 2 of the theorem. However we can obtain 
the desired function by setting 

00 

E(x) = / <j>(xx)exp(-x) dx . 

0 

Clearly E(x) is continuous and satisfies statement 1 in 
A(M). To see that E(x) satisfies statement 2, let x t M and 
t > 0. Then E(xt) £ E(x) holds because <j)(xt) £ <{>(x) holds. 

To rule out E(xt) = E(x), observe that in this case we 

must have <f>(x(t + t)) e <f>(xx) for all x > 0. Thus in 

particular, letting x = 0, t, 2t, ... we get <}>(x) = <}>(x(nt)), 

n = 1, 2, ... But asymptotic stability of M implies that 

for x e A(M) , d(xt,M) + 0 as t + +«. Thus 

<(>(x(nt)) -► 0 as n as <f>(x) is continuous. This shows 

that 4>(x) = 0. But as x t M, we must have <f>(x) > 0 and 

we have a contradiction. Hence E(xt) < E(x) for x £ M 

and t > 0. Thus statements 1 and 2 are satisfied. // 

This theorem says nothing about the size of the region 
of attraction of M. Thus if a function E(x) is known to 
exist in a neighborhood N of M, we need not have either 
N C. A(M) or A(M) C. N. This means that the above theorem 
cannot immediately be stated as a theorem on global 
asymptotic stability. The following example helps to 
illustrate this point. 
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EXAMPLE 6.3 Consider a dynamical system defined in the 
x,y-plane by the differential equations 



# 

x = f (x,y ) , 

y = -y, 



where 



f(x,y) 



2 2 

x is x y > 1 



•5 2 2 2 

2x^y d - x if x y < 1 



The trajectories of this system are sketched below. 





The origin (jj) is asymptotically stable, with the 
{( x ) :x 2y2 < as its region of attraction. Consider 

y 

the function 



set 

now 
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E(x,y) = y 2 + - . 

1 + 

This function satisfies the conditions of Theorem 6.2 in 

the whole plane. To see this we note 
dE _ 3E dx , 3E dy , . , • 

dt “ 97 dt + ay dt which implies E(x,y) = |ff(x,y) - y 

It can be verified that E(x,y) < 0 for all ( x ) ft (°), which 

implies that E satisfies the conditions of the theorem in 

every neighborhood of the origin. But not every neighbor- 
hood of the origin is contained in the region of attraction, 
nor does it contain the region of attraction. 

In view of this example, we wish to establish conditions 
necessary to ensure global asymptotic stability. This is 
done by starting with the following definition. 

DEFINITION 6.4 A continuous real-valued function E(x) 
defined on a set N <1 X will be said to be uniformly unbounded 
on N if given any a > 0 there is a compact set K C N, 

K ^ N, such that E(x) >_ a for x i K. 

The following three theorems will be stated without 
proof since the proofs are very laborious. The proofs may 
be found in Bhatia and Szego, pp. 70 - 73- 

THEOREM 6 . 5 Let M C X be a compact asymptotically stable 
set and let M be_ invariant . Then there exists a continuous 
uniformly unbounded function E(x) defined on A(M) such 
that 
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1. for x e M, E(x) = 0 and for x £ M, E(x) > 0, 

2. for all x e A(M) and t e R, E(xt) = e _t E(x). 



^ • 6 If M C X is any compact asymptotically stable 
set_, there exists a continuous , uniformly unbounded function 
E(x) on A(M) such that 

1. for x e M, E(x) = 0 and for x £ M, E(x) > 0, 

2. for x £ M and t > 0, E(xt) < E(x). 

THEOREM 6 . 7 Let M CL X be_ compact and let there exist a 
continuous uniformly unbounded function E(x) defined on 
an open neighborhood N of M such that 

1. for x e M, E(x) = 0 and for x £ M, E(x) > 0, 

2. for x £ M, t > 0 and x[0,t]C N, E(xt) < E(x). 

Then M is asymptotically stable and N C A(M) . If in 

addition , any condition guaranteeing the invariance of N 
holds , then N = A(M). 

These three theorems are very similar in their conclu- 
sions but we can see that the hypotheses are becoming more 
general. The next results on global asymptotic stability 
follow from these theorems. 

THEOREM 6.8 A compact invariant set MC X is globally 
asymptotically stable if and only if there exists a continu - 
ous uniformly unbounded funct ion E(x) defined on X such that 

1. for x e M, E (x ) = 0 and for x £ M, E(x) > 0, 

2. for all x e X and t e R, E(xt) = e _t E(x). 
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PROOF The sufficiency follows from Theorem 6.7 since X is 
an invariant neighborhood of M. The necessity follows from 
Theorem 6.5. // 

THEOREM 6.9 A compact set M C X is globally asymptotically 
stable if and only if there exists a continuous uniformly 
unbounded function E(x) defined on X such that 

1. for x e M, E(x) = 0 and for x i M, E(x) > 0, 

2. for x l M and t > 0, E(xt) < E(x). 

PROOF The sufficiency follows from Theorem 6.7 and the 
necessity follows from Theorem 6.6. // 

These theorems allow us to determine whether a set is 
asymptotically stable or not without having to solve the 
system. However, if a set is asymptotically stable, it 
may not be a trivial problem to display the continuous 
uniformly unbounded Liapunov function. As in the last 
section, we now focus our attention directly on systems 
of differential equations in order to apply some of the 
above theory. 

Thus we return to the general case of an autonomous 
system 



X = F(X) 

defined in a region of R n with a critical point at Xq 
in fl. Without loss of generality, we may assume that the 
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critical point is the origin 0 of R n for if it is not, a 
simple change of variable X* = X - X Q will translate the 
origin to the point Xq without changing the form of the 
system. 

Again we seek a real-valued function E(X) of class C 1 
defined on some region ft of the origin such that E(X) > 0 
and E(X) = 0 if and only if X = 0. We will refer to E as 
the energy function because of our previous physical 
interpretation. If X = X(X Q ,t) is the solution of 

t 

X = F(X) 

X(0) = X 0 , 

3 E 

the time rate of change of energy, along the trajectory 
defined by this solution is 



n 



VE 



X = VE-P = Z 



3E 



i=i 8x i 



Thus if this expression is negative in ft, a particle moving 
along any trajectory defined by the system will be dissipa- 
ting energy as it intersects the surface E(X) = c. Conse- 
quently a particle which enters the region enclosed by 
such a surface can never gather enough energy to escape, 
and its trajectory must remain in that region. This implies 
that the origin will be asymptotically stable, for the 
given system. 
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These remarks motivate the following definition. 



DEFINITION 6.10 Let A be a region of R n containing the 
origin, and let E = E(X) be a real-valued function of 
class C 1 in A. Then E is said to be positive definite if 

1. for all X in A, E(X) >_ 0 and 

2. the equality E(X) = 0 holds if and only if 
X = 0. 

^ 3 E 

If in addition we have VE*F = E ^ — F. < 0 everywhere 

i=l 3x i 1 “ 

in A, then E is said to be the Liapunov function for the 
autonomous system X = F(X). 

In terms of this definition, we can now form some 
conclusions about the stability of the origin. 

THEOREM 6.11 The origin is a point of stable equilibrium 
for 



X = F(X) , 

F(X q ) = 0 

provided there exists a Liapunov function E for the system . 

PROOF Assume the existence of the Liapunov function for 
the system and consider the set 

S[0,r] = {X e R n : d(0,X) <_ r} . Let r be so chosen that 
S[0,r] lies in A. Since S[0,r] is compact and E is 
continuous and positive definite in A, it follows that E 
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achieves its minimum value on S[0,r]. Let m be this 
minimum value. Clearly m is positive. 

Again since E is continuous, there exists an e > 0 such 



that | | E ( X ) - E ( 0 ) | | < m whenever 
stated, E(X) < m whenever ||X|| < 
Consider the set S(0,e). Let X Q 
be any point in S(0, e ) and let 
X(X 0 ,t) be the unique solution 
of X = F(X) whose trajectory 
passes through Xq at time 
t = 0. 



| | X — 0 | | < e, or simply 




Now we claim that the trajectory 

must remain within S[0,r] for all t > 0. Assume that this 

is not the case. Then | | X( X Q , t x ) | | = r for some t 1 > 0, 

since the trajectory is a connected set. Thus 

E(X(XQ,t^)) >_ m but this implies that E is increasing with 

d x 

time and that > 0. Now = VE-F.' Hence 

9t 9 t 9x. dt 



VE*F > 0. But this contradicts the assumption that 
VE*F < 0 in !1, since this inequality implies that E is 
a non-increasing function of t along the trajectory in 
question. Thus the trajectory remains within S[0,r] for 
all t > 0 and hence the origin is stable since S[0,r] is 
a positively invariant neighborhood of the origin. // 



This argument suggests that if E actually decreases 
along trajectories in a neighborhood of the origin, then 
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the origin would be asymptotically stable. The next theorem 
uses this fact. 

THEOREM 6.12 Lf E is_ a Liapunov function for X = F(X) with 
the property that -VE’F is_ positive definite in ft, then 
the origin is asymptotically stable . 

PROOF We have seen that there exists an e > 0 with the 
property that E is a non-increasing function along any 
trajectory X(X Q) t) with | | X Q | | < e. Since E is positive 
definite in ft, we claim that E approaches a non-negative 
limit Eq along this trajectory as t 00 . To see this, we 
form the sequence {E(X(Xq , t fi ) ) } where t^ -*■ 00 . 

Now E(X(XQ,t k )) assumes values in [0,E(Xq)]. Since 
[0,E(X 0 ) ] is a compact set, every sequence in [0,E(Xq)] 
has a subsequence which converges to a point Eq in [0,E(X Q )]. 
Moreover Eq is non-negative. 

We will be done if we show that the limit is zero when 
-VE • F is positive .definite in ft. Then the fact that E is 
zero only at the origin will imply that the trajectory 
approaches the origin with increasing time. 

On the contrary, suppose Eq > 0. 
continuous there exists a 6 > 0 such 
that | | E(X) - E(0)|| < Eq when 
[ | X | | < 6. Consider the set 
S[0,6]. Now the trajectory 
X(X Q ,t) cannot enter the sphere 
S [ 0 , 6 ] for if it did, then 



Then since E is 




E(X(XQ,t^.)) < Eq for some k > 0. But since E is a non- 



-VE*F is continuous and hence assumes its minimum value on 
the "annular" region 6 5. I l x l I < r where r is as in the 
proof of the previous theorem. This fact is true because 
this "annular" region is compact. Let m be this minimum 
value of -VE*F. Since -VE*F is positive definite, it 
folloitfs that m > 0. Thus 



Thus we have E(X(XQ,t)) E(Xq) ” • 

Since the right hand side approaches -» as t + 00 , this 
implies that E(X(X Q ,t)) < 0 for all t > t Q where 



increasing function, it could never leave S[0,6] and thus 



could not approach E^ as a limit . 

Now since F- is continuous and E is of class in ft , 



3E(X(X Q ,t)) 



< -m for all t > 0 



Now 



E(X(X Q ,t)) - E(X(X q ,0) ) = / 



t 3E(X(X Q ,t)) 



dx 



0 



by the fundamental theorem of the calculus. But 



t 3E(X(X 0 ,x)) 



t 

dx <_ / -mdx = -mt . 

0 



0 



t 




0 



m 



but this contradicts the fact that E is positi 



definite in ft. Thus Eq = 0 and the origin is asymptotically 
stable. // 

Thus far we have presented theorems only on stability 
or asymptotic stability of the origin. However it is also 
important to know under what conditions the origin is 
unstable. Consequently we present the following instability 
theorem known as Liapunov’s First Instability Theorem. 

THEOREM 6.13 Let E be a real - valued function of class C"* - 
in ft and suppose that E(0) = 0, E is_ positive definite and 
E ijj_ able to assume positive values arbitrarily near the 
origin . Then the origin is unstable . 

PROOF Let S[0,r] and S[0,e] be spheres about the origin 
such that e < r and S[0,r] d ft. Since E(X) and the first 
partials of E(X) are continuous in ft, we have that E(X) 
is bounded in S[0,r]. Hence E(X) £ k for all X e S[0,r]. 
Choose any Xq e S[0,e] such that E(Xq) > 0 and let X(XQ,t) 
be the unique solution of X = F(X) whose trajectory passes 
through Xq at time t = 0. 

Now X(X 0 ,t) must eventually cross bdy S[0,r] for if not 
then E(X(X Q ,t)) approaches a limit p in [E(X Q ),k]. Since 
E (X(X Q ,t ) ) k as t -*■ it follows that E(X(X Q ,t)) -*■ 0. 

Now E(X) = VE-F and since E and F are continuous, we have 
that E(X) is continuous. Then since E(X(XQ,t)) + 0, it 
follows that X(X Q , t ) -► 6 since E is positive definite. 

But this is a contradiction since if X(XQ,t) + 0 as t + 
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it must follow that E(X(X Q ,t)) 0 because E is continuous. 
Thus since E(Xq) > 0 we have E(X(XQ,t)) < 0 but this is 
a contradiction. Hence XCX^jt) must eventually cross the 
boundary of S[0,r] and thus the origin is unstable. // 

The important feature of these theorems is that they 
enable us to determine the stability or instability of 
critical points without actually solving the system. This 
is particularly important in cases where the solutions are 
impossible to obtain in closed form or are difficult to 
analyze. Vie can still gain valuable information on their 
stability provided that we can construct a Liapunov function 
for the system. Although there is no general method for 
constructing these functions, we shall see that this is 
fairly easy to do in many cases. To solidify our under- 
standing of this notion, we consider the following example. 

• 2 2 

EXAMPLE The function E(x,y) = x + y is a Liapunov 
function for the system 

. 2 

x = -x + x y, 

2 

y = -y + xy . 

Clearly, E is positive definite and of class C 1 in the 
entire x,y-plane. Moreover, since 

7E = H e l + If e 2 = 2xe l + 2ye 2 and 
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we have 



F(x,y) = ( -x + x y)e 1 + (-y + xy 2 )e 2 



VE-F = -2(x 2 + y 2 ) + 2x 3 y + 2xy 3 = 2(x 2 + y 2 )(xy - 1). 

Hence -VE-F is positive definite in the region xy < 1 , and 
E satisfies all the requirements of Theorem 6.12 in this 
region. Thus the origin is an asymptotically stable 
critical point for this system. 

We shall now describe a method of constructing a 
Liapunov function. To do so we begin by constructing a 
Liapunov function for the constant coefficient linear 
autonomous system 

X = AX 



when all the eigenvalues of A have negative real parts. 

It is noted that we previously determined that the origin 

is asymptotically stable for this system. Thus we will 

not gain anything new from the construction but it will 

allow us to extend some of our earlier results to non-linear 

systems. The construction goes as follows. 

Let e, , .... e denote the standard basis vectors in 
1 n 

R n and, for each i, 1 < i <_ n, let 
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x ± (t) = X jL (e i ,t) = 



denote the solution of this system which satisfies the 
initial condition X i (0) = e^. Then if 

Y = y ] _e 1 + ... + y n e n is any vector in R n , the function 

x(Y,t) = y 1 x 1 (t) + ... + y n x n (t) 

is the solution of the system which satisfies the initial 
condition X(0) = Y. We now set 



E(Y) = / | | X( Y , t ) | | 2 dt, 

0 

and note that if this integral converges, then.E(Y) is 
positive definite. It is clear that E(Y) >_ 0. Also if 
Y = 0, then for all i, 1 £ i <_ n, we have y^ = 0 and hence 
X(0,t) = 0. Thus E(0) = 0. Now suppose E(0) = 0. This 
implies that X(Y,t) = 0. Hence y^ = 0 for all i, 1 <_ i £ n 
since otherwise X^(t) must be zero for some i. But then 
(0) ¥■ e^. Thus we must have that Y = 0. 

The proof that E(Y) is defined for all Y in R n can be 
found in Ostberg, Kreider, Kuller, p. 415. 

Finally we show that E is in fact a Liapunov function 
for this system. We compute the value of E(X(XQ,t)) where 



150 



X(X Q ,t) is the solution of X = AX satisfying X(X Q ,0) = X Q . 
Hence we have 



E(X(X 0 ,t)) 



00 oo 

/ I |X[X ,t),s]| | 2 ds = / | | X(X n ,s + t ) I I 2 ds 

0 0 u 



00 

= / I |x(x n ,u ) || 2 du. 

t U 



Hence 



H (x(x 0' t ) > = It ! I l x( V u) l 1 2 du = -I l x < x o> t) l I 2 

t 

as required since VE-F = 3E/3t along the trajectories of 
X = AX. Moreover, since the real-valued function | | X| | 
is positive definite on R n , E satisfies all the hypotheses 
of Theorem 6.12 and again we have proved that the origin 
is asymptotically stable. 

We now use this method to construct a Liapunov function 
for the following linear system: 

X = -X 

y = -2y. 

This system has the coefficient matrix 

’-1 0 
0 -2 
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with eigenvalues -1, -2. In this case the solution X^t) 
and X£(t) such that X^(0) = e^, X2(0) = are given by 



x 1 (t) = 



-t 




r 1 


e 


, x P (t) = 


0 1 


0 




4sn 

OJ 

1 

<D 

1 



Thus if Y = xe^ + ye^. 



E(Y) = / | | xe -t e, + ye~ 2t e 9 || 2 dt 

0 1 * 



, , 2 -2t . 2 -llt N .. 

= /(xe +ye )dt 



2 2 
x . y 

— + X 



The fact that this function satisfies the hypotheses of 
Theorem 6.12 is easily verified and thus the origin is 
asymptotically stable. 

We now wish to extend these results to non-linear 
systems. For this, we recall the notion of the linear 
approximation which was developed in the last section. 
We recall that in the system 



X = F(X) , 



the function F may be replaced by a function of the form 
JX + G( X) , where J is the n*n Jacobian matrix of F and 
G(X) is small in comparison with JX when ||x|| is small. 
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Now let P be a function of class in a region of 
R n containing the origin, and suppose that F(0) = 0 and 
the real parts of all of the eigenvalues of the Jacobian 
matrix J of F evaluated at the origin are negative. We 
now show that under these conditions the Liapunov function 

E = E(X) for the constant coefficient linear system 

• 

X = JX that has been previously constructed is also a 
Liapunov function for the non-linear system 

X = F(X). 

2 

Now by our earlier results, we know that VE* (JX) £ -||X|| . 

Hence 



VE(X) • F(X) = VE(X) • [ JX + G (X) ] 

= VE( X) ♦ JX + VE(X) • G(X) 

< -I |x| I 2 + | |VE(X)| I | |G(X) I I . 

This last statement is a consequence of the Cauchy-Schwarz 
inequality. Now we also know that ||VE(X)|| £ 2k||x||, 
since VE is a linear transformation on a finite-dimensional 
space. Thus we have 



VE(X) • F(X) £-||X|| 2 +2k||x|| | | G(X) | | . 
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We now use the fact from the definition of a linear 
approximation that lim F i i Lf = 0. But then it 

llxll-o l ' x ll 

G f v \ 

follows that lim r i v ~i r = 0* Hence we can find a 

| | X I I -*0 11x11 

-w-\ 

sphere S[0,6] about the origin in R such that 

| |G(X> | | < ^ | |X| | 

for all X e S[0,<5]. Then within this sphere we have 

VE(X) • P(X) 1 - | |X| | 2 + (2k | |X| | )(^| |X| | ) = -|| l x l I 2 

and indeed E is a Liapunov function for the system 

• 

X = F(X). Moreover, since this inequality also implies 
that -VE*F is positive definite in S[0,6], the hypotheses 
of Theorem 6.12 are satisfied and we have established the 
following result. 

THEOREM 6.14 Let_ F = F(X) be of class C 1 in a region of 
R n containing the origin , and suppose that F(0) = 0. Then 
the origin is asymptotically stable for the system X = F(X) 
whenever all of the eigenvalues of the J acoblan matrix of 
F evaluated at the origin have negative real parts . 

As before, there is a companion result when all of the 
eigenvalues of the Jacobian matrix of F have positive real 
parts. In that case the origin is a point of complete 
instability for X = F(X). 
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These statements verify the comments in Section V 
about the conclusions that may be drawn from a linear 
approximation. Finally we show by means of an example 
that no definite conclusions can be made from a system 
whose linear approximation shows that the origin is a 
stable critical point. 

EXAMPLE Consider the plane autonomous system 

x = y - xf (x,y ) 
y = -x - yf(x,y) 

where f has a convergent power-series expansion in a 

2 

neighborhood of the origin in R and f(0,0) = 0. 

Taking the linear approximation, we have 



X 



JX 



where 




1 

0 



In this case J has eigenvalues ±i so by Theorem 5*8, the 

» 

origin is stable for the system X = JX. 

Now let E(x,y) = x^ + y^. Then 
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VE-F = ( 2 *)- 




y - xf (x,y ) ) 
-x - yf (x,y ) ) 




x 2 + y 2 )f (x,y) . 



Hence for this system the origin is 

1. stable when f(x,y) > 0 in some neighborhood 



3. unstable when f(x,y) < 0 in every neighborhood 



This completes our study of Liapunov functions and 
stability theory in general. At this point, the reader 
should have sufficient background in this theory to be 
able to read and understand further results in this area. 
Specifically recommended for further reading is Chapter 
VIII, C" 1 " -Liapunov Functions for Ordinary Differential 
Equations, in Bhatia and Szego. 



2. asymptotically stable when f(x,y) is positive 
definite in a neighborhood of (q); 
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